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SUMMARY 


i 

This  study  examines  the  thesis  that  a  part  ot  the  human 

physiological  system  can  be  simulated  by  a  suitably  constructed 

« 

mathematical  model.  The  model  employed  derives  from  a  class  of 
mathematical  programming  methods  that  were  originally  developed 
for  representing  complex  military  and  Industrial  activities  and 
have  recently  been  used  to  represent  involved  chemical 

t 

equilibria. 

The  motivation  for  this  research  Is  the  long-range  view 
that  a  successful  mathematical  simulation  of  the  human  system 
or  of  human  subsystems  would  provide  an  Important  tool  for 
biological  Investigations.  A  sufficiently  complex  mathematical 
model — that  lsf  a  model  that  embodies  sufficient  chemical  and 
biological  detail  to  represent  a  whole ,  functioning  human  system 
or  subsystem — could  be  used  to  explore  biological  hypotheses, 
environmental  stress  reactions,  and  Interplay  of  dependent 
subsystems,  and  could  serve  as  a  pedagogical  tool  or  even  as  an 
aid  to  medical  diagnosis. 

Of  course,  the  foregoing  long-range  view  is  an  ultimate 
goal.  For  the  moment,  only  the  techniques,  concepts,  and 
characteristics  of  such  a  mathenatlcal  model  are  being  explored. 

This  paper  presents  the  results  of  a  simulation  of  the 
external  respiratory  function.  Respiration,  and  the  consequent 
gas  exchanges  at  the  lung  surf a  es,  Involves  many  ohemloal 
reactions  and  a  transformation  <  f  venous  blood  into  arterial  blood 

This  activity  was  chosen  as  a  U  tt  case  to  explore  the  feasibility 

\ 

of  constructing  a  mathematical  p  )del  of  a  human  subsystem* 


Basing  their  work  on  the  known  physiological  and  chemical 
aspects  of  this  subsystem,  the  authors  have  constructed  a. 
mathematical  model  that,  when  solved,  yields  values  describing 
the  major  phenomena  of  the  subsystem.  The  values  of  some  thirty 
different  molecular  species  as  determined  by  the  model  are  in 
excellent  agreement  with  observed  values . 

In  the  present  form  of  the  model,  emphasis  was  placed  on 
the  chemistry  and  thermodynamics  of  the  subsystem  rather  than 
on  the  physical  dynamics,  such  as  flow  and  mixing.  These  time- 
phased  phenomena  are  expressed  only  implicitly  at  this  stage. 

In  spite  of  these  and  other  approximations,  this  exercise  appears 
to  demonstrate  that  the  mathematical  art  is  now  capable  of 
representing  such  systems  having  large  and  complex  sets  of 
functionally  interrelated  variables,  and  that  there  is  reason 
to  be  optimistic  about  the  possibilities  of  representing  other 
more  complicated  human  subsystems  and  interconnecting  them. 


FOREWORD 


Traditionally,  the  Doctor  of  Medicine  has  been  trained 
first  in  the  basic  sciences  of  anatomy,  biochemistry, 
physiology,  pathology,  and  pharmacology.  Subsequently,  in  his 
clinical  years  he  has  been  trained  to  integrate  and  relate  the 
knowledge  obtained  in  these  basic  sciences  of  medicine  as  they 
apply  to  the  individual  patient. 

Since  the  turn  of  the  century,  the  accumulation  of 
knowledge  —  not  only  in  the  clinical  fields  but  also  in  the 
basic  sciences  —  has  become  so  vast  that  there  has  neces¬ 
sarily  developed  a  fragmentation  into  many  specialties  and 
Bubspecialtles. 

The  human  body  may  be  considered  in  one  aspect  as  per¬ 
haps  the  moat  complex  chemical  factory  ever  devised.  It  is  a 
dynamic  chemical  factory  in  which  there  are  no  absolute 
values.  The  recent  Introduction  of  isotope "techniques  into 
medicine  has  re-emphasized  the  dynamic  state  of  the  human 

i 

body. 

Recognizing  the  dangers  Inherent  in  oversimplification, 
we  might  nevertheless  say  that  in  medicine  we  are  dealing 
functionally  with  highly  comp]  jx  systems  of  oxidation— reduction 
reactions  and  with  control  mec nanisms  affecting  the  varying 
rates  of  these  reactions.  Of  :he  hundreds  of  parameters 
measured  in  medicine,  we  are,  .n  fact,  measuring  indices  of 
rates  of  conversions  at  a  give  n  time.  The  availability  of 


energy  to  control  these  rate,  may  be  a  critical  factor  In  such 

normal  processes  as  growth  and  aging  —  and  in  the  variations 

from  these  normal  processes,  which  we  may  view  as  diseased 

states  or  pathological  conditions. 

New  approaches  must  be  undertaken  to  assist  in  viewing 

the  human  body  as  a  whole.  These  approaches  must  be  applicable 

to  the  integration  and  evaluation  of  information  and  of  a  very 

large  number  of  parameters  that  have  been  accumulated  but  not 

necessarily  Interrelated  in  the  hundreds  of  special  fields  In 
* 

which  research  has  been  undertaken  in  the  study  of  the  human 
body. 

It  now  appears  feasible  to  use  advanced  techniques  of 
mathematical  programming  and  computers  as  one  means  of  gaining 

greater  Insight  into  the  over-all  complexities  of  the  functioning 

# 

of  the  human  body  in  relation  to  Its  anatomical  structure.  The 
present  trial  study  has  demonstrated  the  feasibility  of  this 
approach  as  applied  to  oxygen  utilization. 

Crawford  P.  Sams,  M.D.,  Research  Physician, 
Institute  of  Engineering  Research 
University  of  California 
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I.  THE  GENERAL  OBJECTIVE 


The  only  way  of  real  advance  In  biology  lies  in 
taking  as  our  starting  point,  not  the  separated 
parts  of  an  organism  and  its  environment,  but  the 
whole  organism  in  its  actual  relation  to  environ¬ 
ment,  and  defining  the  parts  and  activities  in 
this  whole  in  terms  Implying  their  existing  re¬ 
lationships  to  the  other  parts  and  activities.  — 

J»  B.  S»  Haldane 

Haldane's  formulation^  of  the  prerequisite  for 
important  advances  in  biology  was  made  by  that  eminent  British 
biologist  thirty-seven  years  ago.  Until  very  recently,  its 
implementation  had  seemed  well  beyond  the  capabilities  of 
human  accomplishment.  This  is  not  to  imply  that  biologists, 
physiologists,  and  medical  researchers  have  failed  to  develop 
a  knowledge  of  the  micromechanisms  of  the  reactions  occurring 
in  the  separated  parts  of  organisms.  On  the  contrary,  these 
scientists,  working  in  many  countries,  have  amassed  a  re¬ 
markable  amount  of  data  relating  to  reactions  occurring  in 
organisms  and  parts  of  organisms  ranging  from  the  single  cell 
to  complicated  multi functioning  organs  of  the  human  system. 

The  problem  of  interrelating  and  synthesising  these  multitudes 
of  dependent  reactions  into  entire  systems  or  even  into  sub¬ 
systems  representative  of  the  organism  has,  however,  remained 
unsolved. 

The  diagnostician  has  res  >rted  to  intuition  in  attempting 
to  synthesize  in  his  mind  the  results  of  tests  of  isolated 
functions  in  terms  of  the  subs  stem  of  particular  interest. 


Physiologists  have  had  no  suitable  means  for  relating  a  new 
understanding  of  the  functioning  of  a  part  of  an  organism  to 
larger  parts  or  to  the  system  as  a  whole.  Chemists  developing 
new  drugs  have  been  handicapped  in  obtaining  an  appreciation 
of  the  possible  effects  and  side  effects  of  these  drugs  in 
vivo. 

The  coming  nuclear  and  space  age  presents  an  especially 
difficult  challenge  to  the  biological  sciences  and  to  the 
medical  arts.  The  stresses  that  may  be  placed  on  organisms  by 
the  completely  new  and  stringent  environments  associated  with 
these  technologies  are  largely  unknown.  An  understanding  of 
the  effects  of  these  environments  on  the  whole  organism, 
especially  on  the  human  system,  will  be  difficult  and  slow  if 
not  impossible  to  obtain  within  the  present  state  of  the 
biological  arts  and  sciences.  This  is  partly  because  of  our 
lack  of  knowledge  of  individual  biological  phenomena  and 
partly  because  of  the  present  lack  of  a  technique  allowing 
the  integration  of  a  very  large  number  of  environmental 
effects,  which  may  operate  primarily  through  a  part,  on  the 
whole  system.  We  shall  not  try  to  decide  in  advance  which 
lack  is  more  important.  However,  mathematical  techniques  and 
computing  facilities  are  now  becoming  available  for  integrating 
Increasingly  complicated  se*:s  of  components.  It  is  felt  that 

t 

these  techniques  should  be  wohed  to  the  biological  or  mathe¬ 
matical  limit  —  whichever  ome3  first  —  in  the  hope  of 
constructing  mathematical  m  -dels  for  increasingly  large 
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portions  of  the  human  system. 

The  idea  of  building  mathematical  models  o£  biological 
systems  is,  of  course ,  not  new.  Many  Investigators  have 
explored  the  possibilities  in  the  past  and  the  field  is 
attracting  increasing  attention  at  present.  The  names  of 
Lotka,  Rashevsky,  Henderson,  Michaells  among  others  come  to 
mind  as  early  Investigators  who  proposed  mathematical  models 
of  biological  systems  or  who  pioneered  in  the  quantification 
of  parts  of  such  systems.  Present  interest  in  *  \e  subject  is 
attested  by  the  devotion  of  two  complete  issues  of  Reviews  of 
Modem  Physics  (January  and  April  1959)  to  papers  dealing  with 
the  quantification  of  biological  r'ocesses.  Welnrauch  and 
Hetherlngton  have  recently  reviewed  the  increasing  activity  in 
the  application  of  electronic  computers  to  the  study  of 
biological  problems.  (J.A.M.A.  January  17,  1959  PP*  120-125.) 


II.  ON  THE  FEASIBILITY  OF  A  MATHEMATICAL 


MODEL  OF  TKE  HUMAN  SYSTEM 

In  considering  the  forgoing  situation,  the  authors 

believed  that  the  difficult  biological  problem  cf  relating  t'ne 
parts  to  the  whole,  and  the  whole  to  toe  -  uvironroent,  might  be 
well  served  by  employing  some  of  the  newer  mmthemetlcal  pro— 

grarunlng  techniques  that  have  b^en  eppl  e.l  successfully  In 
representing  other  systems  of  comparable  complexity  and  number 
cf  components.  Thus,  linger  ano  nonlinear  programming  rroJ -1* 
h»,vo  been  successfully  constructed  to  scvodule  the  activities 
cf  vs 8 1  industrial  and  military  enterprises. 

The  general  approach  !?  to  view  the  collection  of  organs 
cf  the  body  —  the  heart,  lungs,  liver,  skin,  arteries,  veins, 
.-to.  —  as  art  Interrelated  and  interdependent  physical  set  cf 
subsystems  each  of  which  performs  a  well-defined  function  or 

h.  ;tivity  such  as  th«.-  pumping  activity  of  the  heart.  For  u 

1 

particular  org»n,  the  level  of  activity  met  be  consistent 

4 

with  certain  inputs  foo'T  ot  er  organ?  ana  must  not  exceed  the 
given  capacity  cf  toe  organ  -  for  example,  the  pumping  ac¬ 
tivity  of  the  heart  is  limit'd  by  its  phy  lea’  t^izo  and  the 
f-vailabil  1  ty  of  oxygen  to  he  rt  muscles,  and  also  by  the 
volume  cf  blood  av?  I  Table  fr  w  the  vein?.  The  result  cf  a 
given  level  of  activity  is  c  rtain  output^-  t*  other  organs. 

If  it  wop.  valid  to  e  j  e  that  a  chang*-  in  the  ctivi  ty 
Ir-vel  of  an  organ  results  ir  a  proportion  *  -,hanr.?  in  its  in¬ 
puts  and  outputs,  t*'  n  the  »  -..ultinr  math  *p.-’ tie?’  n-> J-.-l  wou  d 
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be  of  the  pure  linear— programming  type.  It  is  immediately 
evident,  however,  that  gross  functions  such  as  the  respiratory 
function  or  the  liver  function  do  not  remotely  satisfy  the 
condition  of  proportional  inputs  and  outputs.  Indeed,  these 
systems  can  function  under  an  almost  endless  variety  of 
complicated  inputs  producing  a  correspondingly  compli 'ated 
variety  of  outputs.  Nevertheless,  if  we  broaden  our  defini¬ 
tion  of  activity  to  include  the  forming  of  various  chemical 
species  within  some  gross  function,  then  these  more  elementary 

functions  in  most  respects  satisfy  such  an  assumption  of 
0 

proportionality,  the  well-known  chemical  law  of  combining 
proportions.  In  other  respects,  these  elementary  functions 
are  quite  different.  Thu3  the  amounts  of  various  substances 
formed  when  venous  blood  and  air  come  in  contact  through  the 
alveolar  surfaces  of  the  lungs  inu3t  also  satisfy  known  laws 
for  chemical  equilibria.  In  other  cases,  3uch  as  blood 
circulation,  basic  relationships  drawn  from  the  science  of 
hydrodynamics  play  an  Important  role. 

TCius  a  human— system  model,  while  possessing  many 
characteristics  of  those  used  successfully  to  program  human 
group  activitie:  ,  has  other  characteristics  that  result  In 
a  more  complex  structure.  One  aspect  of  the  mathematical 
system  is  clear.  It  must  involve  many  variables,  representing 
levels  of  activity  of  cer“ain  organs  and  rates  of  flow  of 
many  substances,  and  thes  •  must  satisfy  a  large  number  of 


equations  and  inequalities  expressing  the  limitations  of 
capacity  and  availability  from  other  organs  and  sources. 

It  is  the  authors'  belief  that  the  successful  development 
of  mathematical  models  capable  of  handling  such  a  multitude  of 
variables,  and  especially  the  more  recent  mating  of  these 
types  of  models  with  chemical  thermodynamics,  provide  en¬ 
couragement  that  the  complexities  of  physiology  may  be  amenable 
to  mathematical  representation  and  solution. 

Developing  a  complete  and  Integrated  human-system 
representation  is  q  -ite  a  large  order.  In  tne  first  place, 
the  size  of  medical  libraries  al  :>ne  attests  to  the  large 
volume  of  special  cases  that  such  a  model  would  be  required  in 
some  sense  to  interpret.  One  might  nevertheless  argue  that, 
with  luck  and  several  years  of  effort  by  a  large  number  of 
qualified  specialists,  it  might  be  possible  to  handle  the  task 
if  it  were  not  for  one  outstanding  fact.  —  the  great  gap3  in 
our  knowledge  of  the  human  system.  To  cite  one  obvious 

example,  we  know  the  brain  plays  a  major  role  in  the  control 

* 

of  the  physical  pro  -esses  of  the  body,  but  the  exact  func¬ 
tional  relationships  that  define  the  control  can  at  best  only 
be  hypothesized. 

How  then  can  oue  presume  to  build  a  mathematical  model  of 
the  whole  human  system?  The  only  way  of  doing  this  i3,  as  in 
all  science,  through  the  c  whined  use  of  fact  and  hypothesis. 
Thus,  in  areas  where  the  u  iderlying  mechanism  is  unknown,  it 
is  generally  true  that  eitner  X's  or  Y's  hypothesis 


can  be  used  to  fill  the  gap.  Incorporation  of  a 
hypothesis  into  the  model  may  provide  a  useful  way  of  testing 
the  hypothesis  by  comparing  observations  made  under  a  variety 
of  conditions  with  those  predicted  by  the  model  under  similar 
conditions. 

It  is  the  opinion  of  the  authors  that  a  mathematical 
model  of  the  human  system  can  be  built  largely  on  the  func¬ 
tional  relationships  derived  from  basic  physical  and  chemical 
laws.  Such  a  model  is  not  limited  by  the  range  of  empirical 
observational  data.  Rather,  the  mathematical  system  permits 
the  exploration  of  wide  ranges  of  situations  that  may  repre¬ 
sent  extremes  of  environmental  or  pathological  conditions. 

Even  if  at  first  the  results  of  such  mathematical  explorations 
frequently  do  not  correspond  to  reality,  much  useful  knowl¬ 
edge  can  be  obtained  In  determining  the  reasons  for  the  lack 
of  correspondence. 

The  remainder  of  this  paper  describes  the  results,  of 
applying  these  tcchniqueo  to  the  human  respiratory  system* 
Other  human  subsystems  should  Iso  be  amenable  to  representa¬ 
tion  by  the  techniques  described.  Possibly  in  ascending  order 
of  difficulty,  these  might  be  the  following: 

Total  metabolic  sys' am,  inc :uding  digestion, 
nutrition,  and  excr'  tory  functions. 

Synthesis,  storage,  and  conversion  system,  in¬ 
cluding  functions  o  the  liver,  blood-forming 
tissue,  and  muscles 

Hormonal  regulatory  system. 


Central  and  peripheral  nervous  system. 

Mechano-skeletal  response  system. 

It  Is  quite  possible,  of  course,  that  —  as  the  research 
results  unfold  —  thl3  ordering  may  be  hanged  or  different 
groupings  of  components  may  prove  more  expedient. 

It  is  not  known  whether  the  simultaneous  operation  of 
models  of  several  human  subsystems  in  detail  would  be  beyond 
the  capabilities  of  currently  available  computing  machinery. 
If  each  should  turn  out  to  be  no  bigger  than  the  respiratory 
model,  there  would  be  no  problem.  This,  however,  seems  un¬ 
likely.  Nevertheless,  because  of  some  unique  characteristics 
of  the  chemical  syrtems  represented,  and  their  thermodynamics 
and  the  flexibility  of  the  model,  it  appears  to  be  possible 
to  employ  some  tricks  that  may  allow  an  interesting  number  of 
subsystems  to  be  operated  concurrently  with  present  equipment 
Moreover,  in  many  cases  the  additive  qualities  of  the 
chemistry  can  permit  parts  or  all  of  a  subsystem  to  be  aggre¬ 
gated,  and  yet  allow  the  essential  relationships  with  the 
rest  of  the  subsystem  or  with  other  subsystems  to  be  main¬ 
tained.  This  phenomenon  la  something  like  looking  at  the 
human  system  with  a  microscope.  A  small  field  that  is 
being  examined  can  be  seen  in  great  detail,  while  the  rest  of 
the  system,  still  operatin',,  is  involved  only  grossly. 
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III.  A  REVIEW  OP  THE  PHYSIOLOGY  AND  CHEMISTRY 
OF  THE  EXTERNAL  RESPIRATORY  SYSTEM 

To  confirm  the  belief  that  complex  physiological  systems 
can  be  represented  by  mathematical  programming  techniques,  the 
authors  selected  one  important  subsystem  of  the  human  organism — 
the  external  respiratory  system — as  an  exercise.*  This  sub¬ 
system  Includes  the  functions  of  breathing,  the  transfer  of 
gases  to  and  from  the  blood,  the  complex  oxygen  and  carbon- 
dioxide  transport  in  the  arterial  and  venous  blood,  the 
reactions  between  blood  plasma  and  red  cells,  the  reactions 
within  the  red  cells,  the  transfer  of  oxygen  and  carbon  dioxide 
with  the  body  tissues,  and  the  flow  of  blood  through  heart 
action.  An  important  reason  for  selecting  this  subsystem  for 
mathematical  study  is  that  it  has  been  extensively  investigated 
and  a  large  amount  of  quantitative  data  are  available  for 
validating  a  model. 

Figure  1  illustrates  in  a  schematic  fashion  the  gross 
physiologic a]  aspects  of  the  external  respiratory  system  in 
man.  The  ultimate  purpose  of  this  system  Is  to  provide  oxygen 
from  the  atmosphere  to  oxidize  carbonaceous  food,  thereby 
generating  heat  and  supplying  energy  for  performing  work,  and 
to  return  to  the  atmosphere  the  CO.,  (and  some  of  the  water) 
that  is  formed  by  the  .e  reactions  occurring  in  the  cells  of 


#The  internal  resj  ratory  system  relates  to  the  n  *tabolic 
processes  occurring  at  -he  cellular  level.  These  proc  sees 
are  not  represented  »>xj.  icitly  in  the .  present  exercise. 


Nasal  cavity 


Epiglottis 

Trachea 


Lung 


-Bronchus 


Rib  cage 

Pulmonary  artery 
(venous  blood  ) 


Capillary  beds  * 


Systemic  circulation, 
venous 


Pulmonary  vein 
(arterial  blood) 


Heart 

Pleural  cavity 


Diaphragm 


Systemic  circulation 
arterial 


i 


Fig.  I  —  Gross  physiological  aspects  of  the  human 
respiratory  sy  item  (  schemat ic  ) 
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body  tissues.  The  exchange  of  gases  between  the  Mood  and  tne 
external  atrr.ospnere  tal^es  place  in  the  lungs  by  transfer’  across 
the  extremely  thin  walls  of  the  capillaries  and  the  pulmonary 

alveoli . 

The  air  passages  to  and  from  these  exchange  membranes 

include  the  nasal  cavities,  the  pharynx,  larynx,  trachea, 

0 

bronchi,  and  bronchioles.  *nie  lung  tissue  proper  contains 
an  immense  number  of  irregularly  shaped  air  spaces,  called 
alveolar  sacs.  The  terminal  twigs  of  the  bronchial  tree  open 
into  these  alveolar  sacs;  see  Pig.  2.  The  number  of  pulmonary 

p 

alveoli  in  the  lungs  of  man  has  been  estimated  at  7.1  x  10  . 

The  tremendous  surface  area  presented  by  the  alveoli  for  gas 
exchange  results  in  close— to— equilibrium  conditions  between 
arterial  blood  and  sac-gas  composition  except  under  transient 
stress  conditions. 

Ttoe  lungs  play  a  passive  role  in  inspiration  and  expira¬ 
tion.  Their  change  in  volume  during  this  cycle  is  brought 
about  through  change  in  the  capacity  of  the  thoracic  cavity. 

The  change  in  capacity  of  this  cavity  is  caused  by  the  descent 
and  ascent  of  the  dlaphra^n,  and  by  the  movement  of  the  rib 
cage. 

Only  about  one— half  liter,  or  one— tenth  of  the  total  air  ' 
capacity  of  the  lungs.  It  inspired  and  expired  during  ordinary 
quiet  breathing.  This  i  called  the  tidal  air.  Even  after 
the  most  forceful  explra*  ion,  from  1000  to  lfOO  cc  of  air 
remains  in  the  lungs.  Th  space  enclosed  by  the  bronchi, 


Alveolar  sacs 


Fig.  2 — Diagrammatic  representation  of  a  respiratory  unit  (Ref  2) 


bronchioles,  trachea,  larynx,  pharynx,  and  nasal  cavities  Is 
called  the  dead-oir  space.  It  has  a  capacity  of  about  lbO  cc. 
The  air  In  the  air  sacs  Is  termed  the  alveolar  air.  Various 
techniques  have  been  developed  for  measuring  or  estimating  the 
composition  of  this  sac  air,  which  is  in  equilibrium  with 

r 

arterial  blood  and  which  is  different  from  the  external 
atmosphere. 

The  foregoing  characteristics  of  the  lungs,  especially 
the  3mall  tidal  volume  in  respect  to  total  capacity  and  in 
respect  to  the  large  residual  volume,  permit  the  lungs  to 
accomplish  the  remarkable  function  of  damping  out  the  large 
fluctuations  of  sac-eir  composition  that  otherwise,  during 
the  intermittent  inspiration— expiration  cycle,  would  subject 
the  flowing  blood  to  large  surges  of  oxygen  (and  alternately  o 
carbon  dioxide).  To  a  first  approximation,  this  lung  charac¬ 
teristic  permits  the  breathing  mechanisms  to  be  represented 
as  a  large  volume  of  air  of  sac  composition  in  equilibrium 
with  the  blood. 

Actually,  the  breathing  mechanism  and  change  of  air 
composition  from  atmosphere  to  lung  3ac  are  quite  complex  In 
detail.  The  blood  gives  off  carbon  dioxide  (and  usually 
water)  to  tne  lung  sacs  and  pic',  i  up  oxygen.  Thus,  the  air 
leaving  the  sacs  is  of  a  differen*  composition  from  the  in¬ 
coming  air.  On  top  of  this  chang  ,  there  are  other  factors 
such  as  diffusion,  nixing,  the  de  d-cir  spice,  the  alternate 
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inspiration— expiration  cycle,  and  the  blood— flow  rate  that 
all  affect  the  composition  of  3ac  air  In  relation  to  atmos¬ 
pheric  air. 

One's  breathing  rate  can  be  controlled  voluntarily  within 

limits,  but  it  Is  essentially  an  involut  tary  act  controlled 

automatically  by  the  rhythmical  discharge  of  impulses  from  a 

group  of  nerve  cells  in  the  medulla  oblongata  of  the  brain. 

This  respiration  center  is  influenced  largely  by  changes  of 

0 

carbon-dioxide  concentration  In  the  blood.  Changes  in  oxygen 
concentration  can  al30  influence  the  breathing  rate,  but  this 
i3  thought  to  be  an  indirect  influence  through  chemoreceptors 
lying  close  to  the  heart.  These  chemoreceptors  also  appear 

I 

to  be  Influenced  within  limits  by  hydrogen— ion  and  carbon- 
dioxide  concentration,  in  addition  to  affecting  heart-action 
and  blood— flow  rates  and  Indirectly  Influencing  breathing. 

The  gaseous  exchanges  between  the  blood  and  the  tissues 
are  the  reverse  of  those  ta'/ing  place  in  the  lungs.  The 
oxygen  concentration  is  lower  and  the  carbon-dioxide  concen— 
tration  is  higher  jn  the ^tissues  than  in  arterial  blood. 
Therefore,  in  passing  through  th<$*  .systemic  capillaries,  the 
blood  gives  up  oxygen  and  absorbs  carbon  dioxide.  The  compo¬ 
sition  of  the  blood  returning  to  the  lungs  varies,  depending 
on  the  part  of  the  body  it  has  served.  There  is  also  some 
admixture  of  arterial  blood  with  the  venous  as  it  enters  the 
lungs.  As  a  result,  the  b 1 c  >d  entering  the  lungs  —  which  Is 


called  mixed— venous  blood  —  is  higher  in  oxygen  and  lower  in 
carbon  dioxide  than  the  most  degraded  venous  blood,  such  as 
that  coming  from  an  extremity  and  sampled  there. 

The  chemical  phenomena  of  the  external  respiratory  system 
in  which  we  are  especially  interested  may  be  represented  sche¬ 
matically  as  in  Fig.  '>.  Here  the  system  is  considered  or 
consisting  of  three  compartments,  corresponding  to  tie  lung 
sacs,  the  blood  plasma,  and  the  red  blood  cells.  Following 
first  the  path  of  oxygen  from  the  lung  sacs,  we  note  that  in¬ 
asmuch  as  the  concentration  of  oxygen  dissolved  in  mixed- 
venous  plasma  is  lower  than  that  dictated  by  the  partial 
pressure  of  oxygen  in  the  lung  sacs,  oxygen  transfers  from  the 
sacs  to  plasma.  Increasing  its  concent'  ition  there  and  In  red¬ 
cell  solution.  This  increased  oxygen  concentration  causes 
the  reduced  hemoglobin  in  the  red  cells  to  combine  with  oxygen 
until  a  new  equilibrium  is  reached  at  the  higher  oxygen  con¬ 
centration.  Hemoglobin  accounts  for  the  large  percentage  of 
oxygen  carried  by  the  blood;  however,  as  is  also  true  o’ 
carbon  dioxide,  its  solution  in  plasma  and  red  cells  3erve3 
as  the  Important  pathway  in  and  out  of  the  blood. 

Conversely  to  the  case  of  oxyg-u,,  the  concentration  of 
carbon  dioxide  i3  higher  in  mixed— v«  nous  plasma  than  that  dic¬ 
tated  by  the  partial  pressure  of  caroon  dioxide  in  the  lung 
sacs,  so  carbon  dioxide  transfers  fr  m  the  plasma  to  the  -sacs 
until  a  new  equilibrium  is  reached  at  a  lower  carbon-dioxide 
concentration.  But  carbon  dioxide  it  stored  by  the  blood  in 


both  plasma  and  red  cell:*  largely  a3  bicarbonate  ion  (and  In 
association  with  reduced  hemoglobin  in  the  carbamino  form). 

As  the  carbon-dioxide  concentration  is  reduced,  the  bicarbonate 
ion  concentration  is  also  reduced.  This  reaction  in  the  red 
cells  Is  catalyzed  by  the  enzyme  carbonic  anhydrase.  As  the 
bicarbonate  concentration  in  the  red  cells  decreases,  an  im¬ 
portant  reaction  called  the  "chloride  shift"  takes  place 
across  the  cell  membrane;  in  this  reaction,  bicarbonate  ion 
from  the  plasma  transfers  to  the  red  cells,  and  chloride  ion 
from  the  red  cell3  transfers  to  the  plasma.  The  reverse  of 
this  shift  occurs  at  the  tissues  when  carbon— dioxide  concert- 
tratlon  increases  in  the  blood.  Here,  bicarbonate  icn  formed 
in  the  red  cells  shifts  to  the  plasma  and  chloride  ion  shifts 
to  the  red  c^lls.  This  shift  results  in  an  Increase  in  the 
carbon-dioxide  carrying  power  of  the  blood.  The  remarkable 
ion— selective  characteristics  of  the  red— cell  membrane  and  the 
resulting  Glbbs-Donnan  equilibrium  make  this  Important  shift 
possible . 

As  mentioned,  carbon  dioxide  is  also  carried  by  hemo¬ 
globin.  This  hemoglobl n— linked  carbon  dioxide,  while  repre¬ 
senting  only  a  small  part  of  the  tocal  carried  by  the  blood, 
provides  a  significant  part  of  the  -.arbon  dioxide  transferred 
from  tissues  to  lung  sacs. 

The  physiology  and  the  chemist  y  of  the  external  respira¬ 
tory  system  are,  of  course,  much  mo'  «  complex  in  detail  than 
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could  possibly  be  indicated  In  this  brief  description. 

# 

This  Introduction,  however,  should  provide  the  bach-ground 
necessary  for  an  unders tanding  of  the  mathematical  represen¬ 
tation  of  the  more  Important  features  of  thi3  system. 


#An  excellent  detailed  description  of  the  external  res¬ 
piratory  system  is  contained  in  Ref. 


IV.  A  MATHEMATICAL  MODEL  OF  THE  CHEMISTRY 
OF  THE  EXTERNAL  RESPIRATORY  SYSTEM 


DESCRIPTION  OF  THE  MODEL 

From  the  physiological  data  available  and  the  working  of 
the  respiratory  system  as  outlined  above  in  Sec.  Ill,  a 
mathematical  model  was  set  up  to  represent  the  more  important 
of  the  known  interrelated  physiological  functions  and  chemical 
reactions  Involved  In  the  human  respiratory  system.  Actually, 
two  models  were  formulated.  Model  I  uses  a  simplified  concept 
of  the  hemoglobin  molecule,  but  otherwise  is  sufficiently 
detailed  for  the  present  purpose.  Model  II,  described  in 
Appendix  D,  Is  based  on  Linus  Pauling's  theory  of  the  complex 
hemoglobin  molecular  structure. 

Figure  4,  which  Is  a  schematic  Illustration  of  the 
external  respiratory  system,  demonstrates  the  relationship  of 
the  Inputs  and  outputs  of  the  mathematical  model  to  the  actual 
system.  The  a's  and  v's  in  the  figure  refer  to  the  Input 
"elements”  coming  to  the  lungs  in  unit  time  from  mixed— venous 
blood  and  from  the  air.  The  x'3  refer  to  the  resulting 
numerous  molecular  species  in  the  arterial  blood  and  in  the 
air  of  the  lung  sacs  as  determine  i  by  the  solution  of  the 
mathematical  model.  At  the  prese  it  stage,  the  cell  outputs 
are  introduced  into  the  model  In  arras  of  the  composition  of 
the  mixed— venous  blood. 


Cell  output  Cell  input 


Fig  4  —  The  relationship  of  the  inputs  ond  outputs  of 
Model  I  to  the  respiratory  system 
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Model  I  was  constructed  to  provide  an  accounting  of  the 
mass  of  all  of  the  "elements"  Involved.  The  concentration  of 
each  of  the  molecular  species  (dissolved  compounds,  gases, 
Iona)  was  established  In  terms  of  its  mole  fraction.  Since 
the  equilibrium  constants  for  the  many  reactions  were  avail¬ 
able,  It  was  possible  to  Interrelate,  within  the  jnoael,  the 
mutual  effects  of  changes  In  concentration  of  any  of  the 
species.  Further,  from  the  equilibrium  constants  and  mole 
fractions,  it  was  possible,  on  the  assumption  of  Ideal  solu¬ 
tion  behavior  in  the  dilute  solutions  involved,  to  establish 
the  thermodynamics  of  the  system  within  the  model. 

Some  previous  studies  at  F.AND^^  provided  a  technique  for 
the  Introduction  of  these  energetics  of  the  system  into  tne 
computational  model  in  an  'efficient  manner.  This  method  makes 
use  of  the  Gibbs'  chemical-potential  concept  with  the 
minimization  of  total  potertial  energy  of  the  system  defining 
equilibrium,  rather  than  use  of  the  equilibirum.  constants  and 
their  more  conventional  relation  to  free  energy.  The  two 
concepts  are  mathematically  equivalent.  In  the  first  part  of 
Appendix  A,  the  Interested  reader  will  find  a  review  of  the 
basic  chemical  and  thermodynamic  concepts  employed  in  the 
present  development. 

The  model  constructed  to  represent  this  system  is  shown  in 
Table  1.  Briefly,  the  Input  cia  .a  of  the  rr.ocel  consist  of 
given  quantities  of  the  "elemen  s”  of  rr.ixea-venous  blood  and 
air,  shown  on  the  left  In  the  t  ole.  The  .output  data  are  the 


MODEL  I.  THE  EXTERNAL  RESPIRATORY  SYSTEM 
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evaluations  of  the  quantities?  repr?3entl  ng  the  asuo  int;  of 
vartous  molecular  specter.1  found  In  three  compartments  rep¬ 
resenting  lung-aae  air,  blood  plasma,  and  red  ce^s.  The 
numbers  appearing  in  the  various  columns  for  the  Individual 
molecular  species  characterize  the  species  and  are  therefore 
Independent  of  their  quantity.  In  greater  detail,  the  inode: 
consists  of  the  following  portions: 

•1.  An  "Input  Section"  consisting  of  given  amounts 
of  essential  "elements"  "1"  of  the  blood  corresponding  tc 
typical  mixed— venous  blood  values  (measured  In  moles),  plus 
amounts  of  /resh  air  "elements"  "1,"  giving  total  Inputs 
bl  "  ai  *  vi*  (For  simple  chemical  reactions,  the  "elements" 
could  be  defined  as  the  atomic  types  composing  the  molecular 
species.  It  was  more  convenient  for  the  present  application 
to  compose  the  various  possible  molecular  species  appearing 
in  blood  from  e  set  of  "elements"  consisting  of  either  atoms 
or  "groups  of  atoms  In  fl  ed  proportions,"  such  as  H*  or 
OJT.) 

2.  Compartment  I,  label!  "Air  Out,”  with  columns 
corresponding  to  different  sp>  ’l?s  of  lung— cac  gase:  In 

equilibrium  with  the  arterial  plasma  and  red  cells. 

3.  Compartment  II,  labeled  "Arterial  Plasma,"  with 
columns  corresponding  to  cheml  :al  species  appropriate  to  the 
arterial  plasma.  When  the  moc il  is  "solved,"  the  corresponding 
values  are  the  equilibrium  amc  mts  (in  moles)  of  the  species 

in  lung  sacs,  plasma,  and  red  :ells. 


\ 

\ 
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Compartment  III,  labeled  "Red  Celia,"  with  columns 
-  orrcspondirig  to  ch.mlcal  specie,  contained  in  red  cell.  .  A? 
noted  earlier,  Modsl  I  uses  a  simplified  structure.  A  more 
complete  theory  is  contained  In  Model  II;  see  Appendix  D. 

The  mathematical  problem  for  Model  I  Is  t..e  followings 

Por  given  values  (i  *»  1,  2,  . ..,  12),  and  for  the 
detached  coefficient  values  ,  given  in  Table  1,  de¬ 
termine  equilibrium  amounts  / ,  (j  »  1,  2,  ...,  30),  and 

. 1 

minimum  free  energy  z  defined  by 

z  "  f  *i  +  *j  [\i+'n(xy'n)] 

+  III  *J  [CJ+  n*’‘!/TIIl']  * 

where  ^  denotes  the  sum  over  all  species  J  In  compart- 

k 

ment  k  (k  *  I,  II,  or  III),  Cj  Is  proportional  to  the 

free  energy  per  mole  of  the  Jth  species,  x.  ■  2  x,  is 

K  k  J 

the  total  moles  In  compartment  k,  and  In  means  natural 
logarithm,  subject  to  the  constraints 

Xj  ^  0  (j  •  1,  2,  ••*,  30 ) , 

30 

bi  "  alJxJ  ^  *  1*  2* 

Thus,  the  computational  proble  as  presented  In  the  model 

is  to  determine  the  amounts  of  each  of  the  chemical  species 

present  In  the  lung  sacs,  in  arterl-  l  plasma,  and  In  red  cells 
0 

at  equilibrium. 
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In  the  pest,  the  determination  of  equilibrium  compo¬ 
sitions  of  such  large  chemical  systems  Involved  computational 
problems  that  have  b^*en  tedious  and  time— consuming;  they  have 
often  been  "solved”  by  assuming  certain  species  to  be  domlnrnt 
in  the  final  solution. jn  Peoent  years,  however, 
severrl  conceptual  and  computational  developments  plus  the 
availability  of  large  electronic  computing  machines  have  made 
It  feasible  to  attempt  the  solution  of  such  large— ecale 
problems  repeatedly,  and  under  any  desired  variation  of  the 
basic  Inputs  Into  the  system.'  7 

THE  PREg-PnOWT  FUNCTION 

The  method  here  used  for  solving  the  system  Is  described 

in  detail  In  Ref.  5.  It  consists  of  the  direct  minimization  of 
the  nonlinear  free-energy  function  subject  to  certain  linear 
constraints  Instead  of  the  usuil  nonlinear  equillblrum 
equations  that  express  the  nass-actior  laws.  Mathematically, 
the  free  energy  F(X)  of  a  mixture  of  n  chemical  species 
relative  to  some  standard  state  containing  moles  of  the  J-th 
species  can  be  expressed  as 

F(X)  -  £  x  P,  ,  (:) 

J-  J  J 

where  X  •  (x^,  x2*  •  ••#  -n)  1  he  sot  of  mole  nur.be rs , 

represents  the  free  energy  p  r  mole  of  the  J— th  specie' 
at  the  given  temperature  aod  cor  :entretioa#  x./  x;  these  are 
defined  by  *  * 
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_ 

r,  -  RT  (c  ♦  In  -J.  )  , 

x 


n 

Z  x. , 
1  J 


P° 


R*  +  ln  P  • 


(2) 

C) 

(’0 


P  »  pre33ure  ln  atmospheres, 

T  •  absolute  tempers  tv.  re  Ke  1  vi  n  , 
R  »  the  gas  constant, 

P^  ■  free-  energy  per  mo!e  of  the 
J— th  specie.-  ln  Its  standard 

state . 


(Por  species  in  liquid  state,  the  pressure  terra  In  P  defining 
2 j  Is  omitted  for  P's  near  atmospheric;  and  for  pure  species 
in  the  solid  state,  both  thl3  term  and  the  term  ln<  V  x)  are 
dropped. ) 


Where  there  are  several  compartments  In. tend  of  one,  as 
In  the  mode],  it  Is  ne  -t-sary  to  modify  the  free— energy 
equation  so  that  ''oncentr*  tlon  is  defined  »ep.'rately  for  each 
compartment.  The  necessary  change  In  P(X)  1;  ^’.own  in  the 
foregoing  formulation  of  the  mathematical  problem  for  Mode!  I, 
Real  gase:  and  solutions  do  not  conform,  over  the  ntire 
range  of  componli ’ on  and  of  pressure,  to  the  ide* 1  situation 
described  by  these  equations.  To  retain  the  advantages  of 
the  simple  form  of  the  equations  it  has  been  expedient  to 
define  a  thermodynamic  function  c  tiled  activity,  aj,  which  Is 


proportional  to  concentration  In  iole  fraction: 
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The  constant  of  proportionality,  called  the  activity 
coefficient,  approaches  1  a3  x^/  x  tend.:  to  0.  In  this  dis¬ 
cussion,  mole  fraction  x./  x  Is  used  Instead  of  activity  a ,, 

J  v 

with  the  understanding  that  an  approximation  is  thereby  intro¬ 
duced.  The  use  of  mcle  fraction  Instead  of  activity  In  the 
present  model  with  the  very  dilute  ro^utlons  and  low  p^ca.:ure 
Involved  does  not  appear  to  introduce  serious  error. 

The  free-energy  function  is  thu3  defined,  once  the 
numerical  values  of  the  c ,  are  determined  for  the  various 
species.  In  the  case  of  complex  mixtures  of  gases,  extensive 
tables  are  available  (see,  e.g..  Ref.  8)  containing  the  free 
energy  values  F^/(RT)  for  a  large  number  of  chemical  species 
and  over  a  large  temperature  range.  Unfortunately,  these  data 
do  not  Include  the  values  for  many  of  the  chemical  species 
present  in  physiological  systems.  There  are,  however,  equi¬ 
librium  constants  available  in  tee  literature  for  most  of  the 
important  reactions  occurring  in  the  external  reaplratort 
system. 

\ 

Appendix  A  des  jribes  in  d  -tall  the  conversion  of  the 
equilibrium  constants  to  the  c,  valuer  shown  in  Table  1. 

w 

Whenever  jpossiblej,  data  were  drawn  fro r  b<  cl  •  phvsical_c.cemica' 
tables  rather* than  physiological  sources  tc  ensure  their 
applicability  under  various  cor.  1 J  tions  devlat lng  from  those 
of  noimal  individuals. 
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THE  MASS-BALANCE  AND  CHARGE  CONSTRAINTS 

The  determination  of  the  equilibrium  composition  is 
equivalent  to  finding  the  nonnegative  set  of  values  Xj  that 
minimises  F(X)  as  defined  by  (1)  above  subject  tc  several 
sets  of  linear  constraints.  The  first  of  these  are  the  mass- 
balance  constraints, 

n 

«  aij  xj  m  ( 1  ■  1*  2,  •••,  m),  (5) 

<J®1 

where  there  are  m  different  types  of  "elements";  indicates 
the  number  of  units  of  element  1  in  a  molecular  species  J, 
and  b^  is  the  total  number  of  element  weights  of  element  1 
present  in  the  original  inputs.  The  values  of  the  are  the 

entries  in  the  column  below  each  species  in  Table  1. 

Chemically  this  is  equivalent  to  stating  that  the  mass  of  the 
numerous  chemical  species  must  be  in  balance  with  those  of 
the  Inputs.  For  example,  the  first  mass-balance  equation  can 
be  obtained  from  Table  1  by  multiplying  the  unknowns  by 
the  corresponding  entries  in  the  "Op"  row  to  yield 

al  ♦  V1  *  bl  *  X1  ♦  xh  ♦  xi7  ♦  x29  • 

In  addition  tc  the  mass-ba lance  constraints  of  the  model, 
there  is  one  expressed  by  the  "cnarge"  row,  which  indicates 
the  conservation  of  electrical  c!  arge.  It  Is  necessary  to 
Introduce  this  Into  the  model  to  account  for  the  remarkable 
ion-eelective  characteristics  of  the  membrane  separating  the 
plasma  and  red  cells.  The  forum*  stion  of  this  set  of 


constraints  implies  the  Olbt s-Oonnan  equilibrium  that  occurs 
with  such  ion— selective  membranes;  it  also  provides  for  the 
ao-oalled  "chloride  shift"  that  occurs  through  this  membrane. 
This  constraint  may  be  formulated  In  the  present  conte rt  as 

2  x-  -  0,  (6) 

II  J  J 

where  ty  the  constant  appearing  In  the  charge  row,  is  the 
not  number  of  positive  or  negative  charges  present  per  unit 
of  the  charged  epeclee.  The  summation  extends  over  all 

t 

species  In  Oompartment  II. 

In  the  operations  of  the  model  so  far,  the  assumption  ha; 
been  mode  that  there  Is  no  net  charge  difference  in  the  input 
venous  blood,  and  the  charged  input  "elements”  were  carefully 
balanced  In  accordance  with  this  assumption.  If,  as  is  shown 
In  Thble  1,  one  of  the  compartments  with  ions  contains  the 
constraint  that  all  charges  in  the  compartment  sum  to  zero, 
and  If  the  input  charges  sum  tc  zero,  it  is  easy  to  show  that 
the  mass— balance  equations  imply  the  electrical  neutrality  of 
the  other  two  compartments.  Since  the  first  compartment  con¬ 
tains  only  electrically  neutral  species,  (6)  implies  that  all 
compartments  are  electrically  neutral.  Besides  providing  for 
the  representation  of  the  Olbbs^D  mnan  equilibrium  across  the 
cell  membrane  end  the  resulting  c  iloride  shift,  the  charge 
equation  could  be  used  as  a  means  for  forcing  other  compo¬ 
sitional  changes  between  compartm  nte  when  the  charge 
difference  is  not  zero. 


I 
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Some  judgment  entered  into  the  choice  of  the  appropriate 
input  "elements"  and  possible  product  specie'’  present  in 
blood.  As  for  the  'atter,  then.*  wer  *  selected  to  yield  a  fair 
representation  of  the  important  chemical  phenomena  that  have 
been  observed  in  the  extern;!  respiratory  system.  Experience 
has  shown  that  the  model  illustrated  In  Table  1  Is  about  the 
smallest  matrix  in  ter. is  of  rows,  columns,  and  compartments 
required  to  represent  these  more  important  phenomena. 

There  are  also  many  possible  choices  cf  "elements"  to  be 
meet  to  form  the  molecular  species.  The  one  selected  here 
happens  to  be  minimum  In  number  (notice  that  there  are  the 
same  number  of  vectorally  ^dependent  columns  as  th'-r-  are 
re a  trie ting  equations).  The  selection  was  also  made  no  th- t 
ther*3  was  a  physical  Interpretation,  in  terms  of  nonover— 
lapping;  Groups  cf  atoms,  that  permittee  obvious  simple 
number  combinations  to  form  t)  *  various  species. 


whole— 


V.  SOLUTIONS  ON  THE  DIGITAL  COMPUTER 


AND 

SIMULATION  ON  THE  ANALOQ  COMPUTER 

PICT t; L-COMPUTEP  30t  ution’s 

To  solve  actual  numerical  case 3,  the  general  one- 
compartment  (one-phase)  chemlcal-equilib riur  code  developed 
at  RAT'D  for  the  method  described  in  Ref.  5  was  extended  to 
cover  the  multicompartment  case.  The  method  is  iterative  an: 
starts  with  an  arbitrary  composition  of  chemical  species 
Satisfying  the  mass-balance  constraints.  F-y  means  of  the 
assumed  chemical  composition,  a  quadratic  fit  is  made  to 
the  free-energy  function,  which  is  then  minimized.  This 
produces  a  new  chemical  composition  that  also  satisfies  the 
constraints  but  has  a  lower  free  energy.  The  process  is  then 
iterated.  The  method  is  extremely  efficient  and  for  this 
reason  has  come  into  wide  use  for  solving  chemical-equilibrium 
problems . 

Model  I  (which  has  11  mass— balance  equations  and  '0 
chemical  species)  can  be  computed  on  the  IBM  70A  to  four  sig¬ 
nificant  figures  in  approximately  one  minute  with  less  than 
20  Iterations.  The  cost  of  computing  the  output  composition 
for  a  given  set  of  Inputs  from  venous  blood  and  air  composi¬ 
tion  is  about  $r .00.  The  currentl  •  used  digital  <"ode  is 
capable  of  handling  3ystem3  consisting  of  10  times  tne  number 
of  equations  and  variables.  Howev  r,  systems  of  such  size  may 
require  a  pronibitive  number  of  it*  rations  before  an  accurate 


solution  can  be  achieved.  On  the  other  hand,  systems  of  two 
to  three  times  the  present  3lze  are  believed  to  be  economi¬ 
cally  feasible.  So  far  as  the  authors  know,  the  solution 
carried  out  to  validate  the  model  (see  Sec.  VI)  represents 

the  largest  simultaneous,  multiphase  chemical  system  ever 

solved  for  equilibrium. 

Por  Model  I,  the  code  used  only  r ,000  of  the  *2,000  words 
of  high-speed  memory  space  available  (including  space  for 
about  4,000  word3  for  the  code  instructions  and  for  about 
1,000  word*  for  the  problem  data,  namely,  the  Input  values  b^, 
the  energy  constant*  Cj,  etc).  In  addition,  the  machine 
has  magnetic— tape  and  drum  memories.  Thus  it  i3  feasible  to 
consider  having  several  human— subsystem  models  within  the 
machine  at  the  same  time;  these  can  be  cyclically  solved  so 
a  to  feed  input  information  to  each  other.  Por  example,  a 
re  <piratory  model  at  the  tissue  level  could  interact 
with  the  present  external  respiratory  model.  (A  simple  blood— 
flow  simulation  would  be  needed  to  Interconnect  the  two  sub¬ 
systems.)  This  two— subsystem  model  should  give  a  fairly 
complete  simulation  of  the  full  respiratory  function  of  the 
human  body.  One  immediate  application  might  be  to  the  study 
cf  body  reactions  when  breathing  *n  confined  atmospheres 
such  as  are  envisioned  for  space  t; avel. 


v. 

N 
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ANALOQ— COMPUTER  SIMULATION 

In  addition  to  obtaining  accurate  numerical  values  using 
a  digital  computer,  the  authors  explored  the  idea  of  simulating 
the  respiratory  fun  tion  (so  far.  Just  the  hemi  a l—esui librium 
part)  on  the  electronic  circuits  of  an  analog  computer.  For 
this  purpose,  the  RAND  analog  computer  (TOAC)  was  used.  Tills 
is  a  differential  analyzer  originally  made  by  the  Reeves 
Instrument  Company,  but  with  many  modification!) .  The  various 
elements  of  the  machine  Integrators  (which  are  potentiometers) 
were  Interconnected  by  means  of  a  wired  plug  board  that  con¬ 
stitutes  the  Model  I  "program"  for  the  machine  Vo  execute.  A 
series  of  dials  can  be  turned  to  vary  the  inputs  to  the  de¬ 
sired  amounts,  „and  the  circuits  of  the  machine  can  be  scanned 
to  read  the  values  of  any  desired  output  species. 

Computer  solutions  of  Model  I  were  undertaken  with  two 
goals  in  mind.  The  first  war  the  obvious  ne’er.  lty  of  having 
fast  accurate  checks  on  the  mathematical  model;  the  e  ond  was 
the  desirability  of  obtaining,  in  physical  hardware,  a  simu¬ 
lation  the  time-dependent  lystem.  Although  these  two  goals 
aro  not  .^dependent,  the  digital  computer,  with  high  accuracy 
and  reasonably  fast  solution  time,  was  used  to  orapute  the 
equilibrium  conditions  and  to  choc  the  mathematical  model. 

The  analog  computer,  wl tb  less  uce  racy,  but  much  faster  <olu— 
tlcn  time,  can  be  used  eventually,  it  is  believed,  to  simulate 
the  dynamic  aspects  of  the  model.  ITie  Initial  effort  to  this 
end  was  accomplished  in  the  fellow  ng  way. 
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*nie  analog  computer,  as  Its  name  Implies,  Is  singularly 
well  designed  for  the  simulation  of  complex  time-dependent 
systems.  With  the  advent  of  mathematical  programing  methods, 
it  became  possible  to  simulate  complex  systems  for  which  the 
basic  problem  can  be  3tated  in  terms  of  the  extremization  of 
a  given  vector  —  for  example,  the  chemical— equilibrium 
problem  of  Model  I  os  stated  above.  By  means  of  the  methods 
of  mathematical  programming,  an  exact  model  of  the  equations 
for  Model  I  was  constructed  on  the  analog  machine.  An  im¬ 
portant  characteristic  of  this  simulation  is  that  the  entire 
system  is  interconnected  and  hence  interdependent  in  Just  the 
way  the  hypotheses  of  the  physiological  system  are  interde¬ 
pendent.  This  means  that  a  change  In  level  of  concentre tion 
or  activity  of  any  chemical  species  or  component  is 
immediately  reflected  throughout  the  system. 

All  the  important  parameters  of  the  system  were  made 
available  on  the  machine  for  manual  or  automatic  adjustment. 
Once  the  model  had  been  constructed  on  the  machine,  the  input 
or  known  parameters  were  adjusted  by  means  of  potentiometers. 
The  machine  quickly  computed  the  equilibrium  values  for  the 
outputs  of  the  chemical  model.  Ty  t  is,  the  model  on  the 
machine,  with  circulating  currents  came  to  an  equilibrium  in 
which  the  current  level  and  voltage  were  constant,  a 
simulation  of  the  chemical  system.  The  outputs,  or  equilibrium 
values,  could  then  be  read  on  voltr  »ters  or  appropriate  output 
devices  and  t lie  comparison  made  to  he  physiological 
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hypothesis.  For  example,  the  effect  on  the  amount  of  water 
in  the  red  blood  cells  was  studied  by  varying  the  amount  of 
salt  in  the  plasma.  This  was  done  by  continuously  changing 
the  setting  of  the  Na+  and  Cl”  input  dials.  Similarly  the 
effect  of  acidity  on  oxygen  pickup  was  studied  by  varying 
the  H*  input,  etc.  It  is  quite  possible  that  such  a  simu¬ 
lator  with  more  subsystems  would  be  an  excellent  training 
device  to  give  a  medical  student  at  the  controls  a  feel  for 
the  effect  of  this  or  that  change  on  the  functioning  of  the 
hum.' n  system. 

To  make  such  a  system  time -dependent  in  the  most 
elementary  way,  the  Input  parameters  (potentiometers)  —  which, 
for  example,  set  the  CO^  concentration  in  the  lung  alveolar 
sacs  —  were  made  to  vary  in  a  sinusoidal  fashion  with  a 
frequency  comparable  to  the  breathing  rate.  The  entire  system 
oscillated  according  to  the  interdependent  relations  at  V o 
same  rate.  These  oscillations  may  not  correspond  to  nature, 
of  course,  because  a  more  complete  representation  may  require, 
for  example,  the  simulation  of  the  time  lag  encountered  acros. 
physiological  boundaries,  the  interconnection  of  lmllar  ; inf¬ 
lations  of  other  organs,  and  a  simulation  of  the  respiratory 
function  at  cellular  level. 

The  simulation  of  Model  I  requ  red  appro  imately  0 
computing  amplifiers  and  0  potentl  meters.  The  small  analog 
computer  at  The  RAND  Corporation  ca  be  made  to  simulate  a 
system  about  twice  th”  size  of  th  •  j  resent  Model  I 


e.g.  , 
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the  full  respiratory  function.  Por  larger  system?,  tnere  are 
machines  available  outside  RAND.  It  is  alro  possible  that, 
as  subsystems  become  fully  validated,  special-purpose  analog 
omputers  could  be  built  for  the  solo  purpose  of  simulating 
the  particular  subsystem:. 


; 


VI.  VALIDATION  OP  MODEL  I 


Where  possible,  the  respiratory  model  used  equilibrium 
constants  derived  from  physi;al  'hemical  data  rather  than 
those  derived  from  laboratory  mea  .urements  on  the  biological 
medium  itself.  The  rationale  for  this  choice  is  that  data  of 
the  second  type  may  reflect  the  presence  of  other  substances 
and  reactions  and  may  be  valid  only  for  the  normal  or  measured 
conditions,  on  the  other  hand,  the  use  of  basi  c  data  of  the 
first  type  hould  allow  explorations  with  the  model  of  con¬ 
ditions  that  are  far  from  normal  in  respect  to  temperature, 
pressure,  blood  and  air  composition,  etc.,  because  the  more 
basic  equilibrium  constants  can  be  expressed  functionally  in 
terms  of  these  conditions.  The  important  first  question  to  be 
asked  is  how  well  a  model  constructed  in  this  manner  conforms 
to  the  actual  physiological  phenomena. 

To  validate  the  model,  the  problem  chosen  was  to  find  the 
equilibrium  amounts  of  each  of  the  species  for  the  inputs 
shown  in  Table  2  for  an  average  resting  adult  male  (at 
atmosphere  and  °C).  Tills  was  accomplished  by  minimizing  the 
sum  of  the  Individual  partial-moia 1  free -energy  functions  for 
the  species,  subje  t  to  the  linear  constraints  for  the  con¬ 
servation  of  mass  and  of  charge,  using  the  electronic  computer 

code  described  in  Sec.  V.  Tab^e  't>  lists  the  calcu  ated  value 

(9) 

from  the  model  and  the  observed  values  reported'  '  for  average 
resting  individuals.  It  can  be  seen  that  the  correspondence 
of  the  calculated  and  observed  values  is  quite  close. 
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I! 


From  Air  (moles  ! From  Mixed— venous 
per  0.0^1  moles  -lood  (mole6  per 
Cl err.ent"  | _ a lr) _ liter) _ _ To tal 


c.: 

,al 

-  7 ' 

.7  /  c3 

x  10“' 

V1 

r; 

6 . f  3  x 

»-» 

0 

1 

V4 

bl 

rt 

1 .  ’  51  7  2  v  10“ 

■:cv 

a. 

C2 

0  t  p  7  v 

l0-< 

v0 

2.27  x 

_o 

10 

ifc 

a 

.  270“  2"'  x  10 

? 

c 

V 

*  r 

C 

,  a 

✓ 

= 

2.441  r 

x  10“  2 

V5 

= 

4.37  x 

10“* 

b:> 

= 

2.4^46  x  10“: 

c 

0 

v4 

a 

46.7 

b 

s 

■•4.70 

CH~ 

1  a_ 
n 

= 

0 

1 

v5 

zs 

46 .71973 

bc 

3 

4,  .716” 5 

:i~ 

af 

a.' 

c 

1 

V 

* 

0.0-n, 

bo 

0.0614 

Na  + 

'  'ly  3  0 

|  Vy  «  C.08092 

by  «=  0.08092 

K4 

1 

ae  =  0 

Vc  =  C.050 

ib  «=  O.O^C 

H’  1"* 

a9  "  0 

v~  *  ;  .09  x  ic 

1  y 

jb.-  *  9.09  x  10“' 

1  ✓ 

v.?r~h 

P 

;,10  ’  0 

■v^q  »  c .  8  0  x  10 

|b^Q  »  .  0  x  10  J 

nr  “b 
r  1 

au  -  0 

v11  -  1.19  X  10“? 

!bal  -  1.19  x  10-2 

, 

Charge  z, 

a,  ^  =  0 

1  £ 

vi2  ■  0 

b12  "  0 

o 


tr.ey  ror ne sent 


These  dato  have  been  taker,  from  -ef. 
averages  for  resting  adult  males. 

c The  symbol  nl  denotes  1  ^  reduced  hemoglobin  molecule, 
mol .  wt.  1  , SCO;  symbol?  r'r  ana  I  represent  miscellaneous 


protean  .  t 


j.lasn. 


V  1 


cc 


:cll 


re  r 


»  4-  4 


COMPARISON  OP  VALUES  OBTAINED  PROM  MODEL  I  AND  OBSERVED  BLOOD  VALUE 


Value  from  Model  Observed  Value® 

Species _ Quantity _ (moles/vol.  produced)  (moles/vol.  existing) 

_  Compartment  III.  Red  Cells _ _ 
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Some  Interesting  features  of  thi3  equilibrium  model  may 
be  observed  by  following  some  of  the  reactions  and  solutions 
in  Tables  1,  ?,and  ; .  For  example,  it  is  important  to  ob¬ 
serve  the  respiration  quotient  (RQ) ,  the  moles  of  C02  out 
divided  by  the  moles  of  0o  in.  The  moles  of  C0o  out  may  be 
obtained  by  subtracting  the  a2  value  for  C02  input  (Table  ?) 
from  0.0326879  times  the  x2/  x  value  for  CO^  output  (Compartment 
I,  Table  3).  If  a  similar  subtraction  Is- carried  out  for  02, 


the  results  give 


°-°Pm876  .  0-56> 

0.008067* 


a  value  well  within  the  range  for  resting  individuals. 
Similarly,  the  per  cent  of  hemoglobin  saturation,  as  de¬ 
termined  by  the  model,  is  the  reasonable-appearing  value 

100  x29  0.878  _  c 

*2Q~+~%  m  ffnwgo?  96*°  * 

Inasmuch  as  the  carbaalno  reaction  involving  the 
combination  of  hemoglobin  and  C02  was  not  represented  In 
Model  I,  It  Was  necessary  to  compensate  for  this  lack  of 

a  C02  sink  in  arterial  blood  by  reducing  the  total  venous  C02 
input  by  a  comparable  amount.  Consequently  the  v2  input  value 
indicated  is  1  •  10“  mole;  less  than  the  average  true  value 
for  total  C0?  ir  venous  blood  to  compensate  for  the  lack  of 
•arbamirie  repr  1  ntition.  Even  though  the  carbuslnc  reaction 
involve?  a  sma’’  T  part  (approximately  '  per  cent)  of  the  C02 
carried  by  the  ’ood.  It  accounts  for  a  much  hi  ;her  portion 


(approximately  yC  per  cent)  of  the  C02  transferred  by  the 
blood  during  the  respiratory  cycle.  A  model  of  the  next 
degree  of  sophistl nation  must  therefore  account  for  this  im¬ 
portant  reaction.*  The  representation  of  the  amounts  of  the 
other  gaseous  constituents  in  Compartment  I  — -  namely,  and 
H^O  —  appears  to  be  nicely  in  line  with  the  inputs  and  with 
reported  values. 

In  the  cases  of  columns  xr>,  x^,  x^,  x^v,  x^g,  x^,  which 
show  the  amounts  of  0-,,  C0?,  and  N?  dissolved  in  plasma  and 
red  cells,  respe  ctlvecy,  it  will  be  noted  that  there  are 
some  rather  wide  variations  between  the  values  calculated 
by  the  machine  and  the  observed  values.  This  is  most 
apparent  for  02  and  N2  dissolved  in  red  cell3.  There  is 
a  possibility  that  the  model  results  may  be  more  nearly 
correct  than  the  observed  values  in  the-e  cases.  Such  very 
small  quantities  of  dissolved  gases  are  difficult  to  de¬ 
termine  in  the  laboratory.  Then,  too,  it  is  possible  that 
small  amounts  of  hemoglobin  remained  in  the  analyti cal 
samples,  which  combined  with  oxygen  and  thereby  made  the 
solubility  appear  to  be  greater  than  it  is  in  the  complete 
absence  of  hemoglobin.  Although  the  number  of  moles  of  these 
dissolved  gases  is  different  for  plasms  and  red  cells,  their 
mole  fractions  are  identical  in  both  compartment.?  in  the 
machine  solution.  ®iis  is  °o  for  all  specie:-,  appearing  ev.c  pt 


*For  an  e>:ci  llent  description  of  the  phytic  ogical  role 
of  carbamino  co&j  rands  in  blood,  see  Ref.  T,  pp.  B'fSH, 
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for  those  charged  species  that  are  Influenced  by  the  Gibbs— 

Donnan  equilibrium.  It  should  be  ncted  that  the  solubility 

coefficients  used  in  the  model  are  for  water,  not  for  plasma 

0 

and  red  cells. 

It  i 8  also  interesting  to  note  that  the  model  creates  < 

water  in  the  three  compartments  from  the  inputs  H+  and  0H~ 
operating  through  the  reaction 

H*  +  0H“  ^  H^O  . 

The  amount  of  water  formed  is  Influenced  slightly  by  the 
appropriate  demands  and  supplies  of  H+  and  OH-  by  the  com¬ 
peting  reactions. 

The  greatest  difficulties  in  connection  with  the  valida¬ 
tion  of  the  model  arose  in  connection  with  the  representation 
of  acidity  and  concentrations  of  H+  and  Otf~,  columns  x^,  x^, 
x?q,  and  x^ ,  in  plasma  and  red  cell  3.  Acidity  of  blood 
samples  is  usually  determined  by  the  use  of  a  pH  meter  and  a 
glass  oi  olnilar  electrode.  It  is  common  tc  **0*  the  ex¬ 
pression 

pH  5  -1  [H4]  or  pH  =  -log]c  . 

However,  in  recent  years  conclusive  arguments  have  been  pre¬ 
sented,  based  on  both  theoretical  and  practical  electrolytic¬ 
cell  consider? tl onr ,  that  pH  values  are  not  measures  of 
hydrogen— ion  concentration  or  activity.  Rather,  pH  determina¬ 
tions  are  a  reliable  measure  of  the  "aridity"  of  the  unknown 
solution  in  relation  to  <:n  empirical  standard  solution  having 


the  same  general  characteristics,  including  temperature. 
Conventional  pH  measurements  are  of  use  primarily  for  control 
and  gross  comparison  purposes,  but  are  not  Suitable  for 
thermodynamic  calculations  or  mass-action  law  use.  There  are 
not  even  really  satisfactory  methods  for  approximating  [H4] 
or  [OlT]  from  pH  values  or  vice  versa. 

The  problem  this  presents  In  the  validation  of  the 
operation  of  the  model  Is  to  determine  If  the  model  calcu¬ 
lates  [H4]  and  ton  accurately  when  the  only  laboratory 
measure  of  these  concentrations  for  comparison  are  pH  values 
for  plasma  end  red  cells.  One  approximation  suggested  in  the 
literature^ is  that 


pH  measured  —  mpH  »  -log-jQ  [H*3  , 

when  rn  Is  determined  from 


*w 


[H,0] 


where  is  established  by  other  means,  such  as  conductivity 

measurements.  It  will  be  noted  that  if  the  reported  pH 
values  for  plasma  and  red  cells  art  taxen  equal  to  — log^  [H4]. 
and  determined  by  subtraction  from  log  Kw 

(which  Is  —13.35  at  y  °C) ,  then  the  resulting  concentration 
product  does  not  equal  K^j  this  indicate;  that.  In  blood, 

pH  measured  ft  — log1^  [H4]. 

However,  if  one  assume:'  that  the  acid— bare  reactions  finally 
involve  only  the  vfatcr  '.o\  talned  in  the  plasma  cno  red-^ell 


-  "Ic' 

i.»vj 

f- 
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material,  and  If  one  redunes  the  given  pH  valuer  (and 
cal elated  tom  values)  to  correspond  to  those  for  which  . 
water  only  Is  present,  then  the  values  so  obtained  for  real 
blood  satisfy  the  several  mass-action  relationships  quite 
nicely  and  compare  well  with  the  values  calculated  oy  the 
model  and  listed  in  Table  2.  This  correction  wss  applied  to 
the  values  of  cm  and  tom  listed  as  observed  ‘values  in 
Table  2  and  as  approximated  from  the  reported  pH  measurement: . 

The  treatment  of  species  x^,  x.,  -,  xOT,  and  x.,q,  sodium 
icn  and  proteins  in  plasma  and  potassium  ion  and  proteins  1\ 
red  cells,  is  described  in  Appendix  B.  Species  and  x?  , 

HCOT  in  plasma  and  in  red  cells,  is,  of  course,  a  crucial  ion 
in  the  blood  both  for  C02  transport  and  for  buffering.  It  is 
created  by  the  model  in  Compartments  II  and  III  by  the  reaction 

C0?  4  OH"  HCO" 

in  the  face  of  the  demands  and  supplies  of  these  Inputs  from 
competing  reactions.  Bicarbonate  ion  is  said  to  hydrolyze 
(see  p.  219  of  Ref.  3.1) : 

HCO"  ♦  H20  HoC0t  4  OH" ; 

and  it  can  ioni/e  further: 

HC07  CO*  4  H*  . 

j  > 

The  combination  of  these  two  reactions  gives 
# 

2HC0T  4  H20  J^CC  .  4  CO"  4  H+  4  OH" . 


4 


*\  - 
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Species  xl4,  x^,  x2C, 


x 2 ry  — -  nemely  H-,CO-> 


and  C0-.  In  plasma 


and  red  cells  —  with  their  appropriate  thermody actinia  func¬ 


tions,  were  Introduced  In  the  model  to  explore  their  effects 

on  concentrations  of  HCOT#  H4,  and  OK".  The  effects  of  HCC" 

Z>  ; 

hydrolysis  and  second— stage  ionization  on  [HCOp  and  acidity 
appear  to  be  minor. 


The  concentrations  of  species  H2C0  and  CO,  are  not  re—  ' 
ported  In  the  literature.  These  species  were  retained  in 
the  model  tc  Illustrate  the  possibilities  of  using  the  mod  1 
for  explorations  of  this  sort.  In  addition,  it  has  heon 
noted  that  the  model  is  extremely  and  correctly  sensitive  and 
responsive  to  changes  in  the  inputs  or  in  the  environment 
(e.g.,  temperature  and  pressure). 
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VII.  POTENTIALITIES  AND  LIMITATIONS  OF  MODEL  I 

The  authors  hasten  to  admit  that  the  model  described  is 
deficient  in  many  respects,  for  it  represents  but  a  first 
attempt  to  determine  the  feasibility  of  programming  a 
complicated  physiological  system  in  this  manner.  Because  the 
interest  of  this  particular  study  is  centered  on  the  equili¬ 
brium  phenomena  of  the  system's  chemical  aspects  rather  than  on 
the  physical  flow,  mixing,  etc.,  several  expedients  have  been 
used  to  avoid  incorporatir -  a  representation  of  the  physical 
and  transient  aspects  of  the  entire  external  respiratory 
system  in  the  model.  For  example,  the  mechanism  of  breathing 
with  its  complex  intermittent  flow,  mixing,  and  diffusion 
effects  has  been  avoided  for  this  initial  study  in  several 
different  ways: 

1.  By  m  lng  the  set  of  input  "elements,"  a^,  ...,  a^, 
representative  of  a  large  excess  volume  of  air  having  the 
observed  composition  of  the  gases  of  the  pulmonary  alveoli  of 
a  normal  resting  male  at  sea  level.  The  assumption  introduced 
by  this  procedure  Is  that  the  concentration  of  gases  in  the 
alveoli  of  resting  individuals  does  not  change  significantly 
during  the  breathing  cycle. 

2.  By  relating  the  required  incremental  input  of  "fresh" 
air  of  atmospheric  composition  to  the  increments  of  gaseous 
"elements"  needed  to  convert  a  unit  volume  of  mixed— venous 
blood  to  arterial  blood.  The  "fresh  air"  inputs,  a,,  . ..,  aR, 
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per  liter  of  blood  may  be  estimated  on  this  basis  when  this 
Information  Is  not  provided  by  the  source  data. 

3.  By  programming  the  computer  to  sample  the  concentra¬ 
tions  of  02  and  C02  in  either  the  arterial  blood  or  sac  air 
compartments  and  calculate  by  a  series  of  iterations  the 
amounts  of  "fresh"  air  of  given  composition  required  to  make 
these  lung  or  blood  values  correspond  to  previously  established 
standards. 

All  three  of  these  techniques  have  been  used  and  are 

satisfactory  within  their  limitations.  Method  (2)  above  w as 

used  in  obtaining  the  values  shown  in  Tables  2  and  3.  These 

devices  avoid  the  necessity,  at  this  early  stage,  of  combining 

a  physical  model  of  breathing  with  the  present  chemical  model, 

which  may  require  the  introduction  of  a  number  of  concepts  in 

« 

addition  to  those  needed  for  the  chemical  model. 

Another  assumption  implicit  in  the  model  is  that  the 
arterial  blood  leaving  the  lungs  comes  to  equilibrium,  through 
the  alveolar  and  capillary  walls,  with  the  gases  in  the  alveoli. 
As  noted  earlier,  realization  of  this  equilibrium  situation 
appears  to  be  closely  approximated  by  resting  individuals. 

Hence,  for  the  one-cycle  case  used  to  validate  the  model 
(approximately  0.55  liter  of  arterial  plasma  plus  approximately 
0.45  liter  of  arterial  red  cells  in  equilibrium  with  each  other 
and  with  approximately  0.033  moles  of  alveolar  gases,  all 

# 

One  mathematical  representation  of  the  breathing 
mechanism  can  be  found  in  ;  ?f.  12. 
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subject  to  the  mass-balance  constraint  provided  by  the  "elements'* 
of  one  liter  of  mixed-venous  blood  pljs  0.021  moles  of  alveolar 
gases),  the  assumption  Is  made  that  equilibrium  Is  reached 
between  the  solutions  within  the  red  cells  and  the  plasma.  A 
comparison  of  the  results  obtained  from  the  model  with  those 
obtained  from  laboratory  data  on  blood  samples  shows  that  this 
assumption  does  not  appear  to  introduce  any  serious  error  for 
resting  Individuals. 

The  multlstaged  Ionization  of  hemoglobin  and  the  other 
proteln-llke  blood  substances  is  not  represented,  and  the 
important  carbamino  reaction  is  only  implicitly  represented. 

In  Model  I.  In  Appendix  D,  an  expanded  system  (Model  II)  Is 
shown,  demonstrating  how  these  factors  may  be  Incorporated. 

At  the  present  stage,  the  model  is  one  of  Instantaneous 
reaction,  with  time  not  explicitly  introduced  as  an  independent 
variable.  This  means,  of  course,  that  the  present  representa¬ 
tion  can  deal  only  with  features  of  the  respiratory  system 
that  are  capable  of  treatment  by  reversible  thermodynamics; 
however,  this  limitation  still  allows  many  matters  of 
importance  and  interest  to  be  explored.  In  general,  if  the 
number  of  relatively  slow  chemical  reactions  entering  the  model 
Is  not  large,  it  may  be  possible  to  develop  a  dynamic  model  by 
introducing  some  simple  delays  in  the  electrical  analog  of  the 
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current  model.* 

It  will  also  be  necessary  to  develop  and  Introduce  the 
feed-back  loops  by  which  changes  In  C02  concentration,  and 
possibly  hydrogen-ion  concentration,  influence  the  ventilation 
rate  and  the  rate  of  blood  flow  through  actions  on  the 
cardlovasoular  system.  Consideration  must  also  be  given  to 
the  boundary  conditions  essential  for  the  continued  functioning 
of  the  subsystem.  One  important  limit  is  the  maximum  osmotic 
load  tolerated  by  cells.  Fortunately,  the  model  is  nicely 

a 

suited  to  the  incorporation  of  this  limit. 

In  spite  of  these  difficulties,  it  should  be  possible, 
through  the  use  of  models  of  this  type,  to  demonstrate  the 
detailed  interactions  of  the  human  respiratory  system  a3  the 
external  environment  changes  and  as  the  inputs  from  related 

subsystems  change  (e.g.,  in  metabolism).  Changes  in  total 

< 

respiratory  gas  pressure  or  relative  mixtures  will  be  reflected 
in  the  model  by  the  appropriate  changes  in  oxyhemoglobin 
concentration,  alkali  reserve,  blood-flow  rate,  breathing  rate, 
plasma  composition,  etc.  Because  the  energetics  are  an 
integral  part  of  the  model,  changes  in  work  load,  temperature. 


In  highly  dynamic  situations,  such  as  when  the  respira¬ 
tory  system  is  stressed  by  exercise,  there  is  good  reason  to 
believe  that  the  controlling  reactions  will  be  those  Involving 
membrane  transfers.  There  are  three  such  transfers  in  the 
external  respiratory  system,  involving  (1)  the  alveolar  tissue 
in  the  lungs,  (2)  the  red-cell  membrane,  and  (5)  the  body-cell 
membrane.  The  introduction  of  the  appropriate  kinetic  equations 
for  any  or  all  of  these  three  barriers  may  suffice  to  explain 
all  the  chemical  time  effects.  The  remaining  reactions  occur  as 
homogeneous  reactions  in  dilute  solution  or  in  the  gas  phase  and 
are  thus  generally  very  fast  compared  to  membrane -diffusion 
processes. 
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and  pressure  Imposed  on  the  human  system  should  be  reflected 
by  the  appropriate  and  detailed  response  within  the  modeled 
respiratory  system.  Changes  In  fields,  either  gravitational, 
electromagnetic,  or  electrostatic,  can  also  be  introduced 
through  the  thermodynamics.  Abnormal  changes  in  any  of  the 
constituents  of  the  respiratory  system,  as  brought  about  by 
pathological  conditions  or  by  unusual  environments,  should 
elicit  the  appropriate  response  In  the  other  constituents, 
including  the  energetic  aspects.  Thus  an  external  respiratory 
submodel  aftfer  further  development  should  prove  most  useful 
for  diagnosis  and  for  predicting  the  possible  effects  of 
unusual  stresses  on  the  respiratory  system  without  endangering 
human  life. 


P-1811 

53 


Api*? no lx  A  # 

A  REVIEW  OF  CERTAIN  CHEMICAL  THERMODYNAMIC  CONCEPTS 


For  the  benefit  of  readers  who  are  not  familiar  with 
chemistry,  we  shall  review  some  of  the  basic  concepts  used 
in  our  study,  especially  the  relation  of  the  mac-' -action 
law  to  minimum  free  energy.  Suppose  a  compartment  contains 
various  species  such  as  HgO,  0JT",  H+,  C02,  HCO,,  when 
CO-  is  dissolved  in  the  solution.  The  amounts  of 

C  • 

these  various  constituents  depend  on  tint-  total  amounts  of 
water  and  carbon  dioxide  originally  ’’dumped"  into  the  compart¬ 
ment,  but  their  relative  concentrations  to  each  other  satisfy 
certain  conditions  that  are  Independent  of  the  input  amounts . 
In  general,  the  concentrations  {or-  activities )  cf  ti  e  pecies 
at  equilibrium  have  the  property  that,  independent  of  the  in¬ 
put  quantities,  one  or  more  ratios  cf  products  of  concentra¬ 
tion  among  them  have  fixed  values.  These  nonlinear  conditions 
are  called  mas3-ectlon  relationships. 

For  example,  if  we  let  x  denote  the  total  number  of 
molecules  in  the  mixture,  so  that 

*  >:K^0  *  X0H  +  'F  +  :C0  +  X TIC 07  1 

where 


v 


th  -  number  of  ionla  •m!  vojt-r  molecule  , 
th  *  numb  r  of  1  ->ns  of  0T  ,  <  . , 


^nd  define 
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H.G 


N 


0} r 


%o 


X0H“ 


,  etc.. 


to  be  the  concentrations  of  these;  species  in  the  mixture,  then 
the  mass-action  law  states  that 


N 


Ho0 


-  K, 


£wr)  (V)  1  ’ 

where  has  a  fixed  value  (at  constant  temperature  and 
pressure),  called  the  equilibrium  constant  for  the  reaction 

±  Ho0  . 


H+  +  Off 


Even  if  no  CO^  were  present,  or  if  a  salt  were  added  tc  the 
mixture,  this  same  relationship  would  hold  oespite  the  fact 
that  the  individual  concentrations  forming  t’.e  above  ratio* 
could  be  changing  drastically.  Another  ma is-ar^ion  relation¬ 
ship  that  holds  among  the  concentrations,  indep  ndent  of  the 


“  Kp  » 


where  K,  is  the  equilibrium  constant  for  the  reaction 

CO.  +  H,  C  ?  hc  or  +  H+  . 
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“Die  conventional  method  for  determining  the  equilibrium 
composition  is  to  find  the  value?  of  the  five  unknowns  :<^  that 
satisfy  five  equations  —  two  of  which  are  nonlinear  and  are 
the  two  mass-action  conditions  given  above,  while  th*.  remain¬ 
ing  three  are  linear  and  express  the  mass  balance  in  relation 
to  the  original  Inputs,  for  example. 


xHoO  +  X0K  +  xHCO- 


xH20  +  XH+ 


*CO„  +  xHCO, 

c. 


input  CH  , 

input  H+, 

inpUt  OOry  . 


As  noted  earlier  in  this  research  memorandum,  we  have  found 
it  more  convenient  to  bypass  the  mass-action  relationships -by 
finding  the  equilibrium  composition  through  direct  minimiza¬ 
tion  of  the  free-«nergy  function#  subject  only  to  the  mass- 
balance  conditions. 

In  general,  for  a  chemical  reaction  represented  by 

P  ^  +  A^,  P  o  4-  .  .  .  A^  P^  +  A,}  P|  +  .  .  •,  ( 7 ) 


the  corresponding  mass-action  relationship  is 


species , 


P,  err*  the  reactant  .speller,;  th;  P,  ere 

J  •' 

i 

the  A.  end  Af  are  th  corresponding  nuir.be 

iJ  J 


the  product 
ri  of  moles 


of  ?ach  -.Dccier; ,  K  1  t\  *  *oui*  Lbriu..  -onr tant,  jnc  the 
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i 


*  T 

l*  .  f 

J 

nd  ar-.~  t’.' 

n:  n  centr-  tioac 

•f  th 

'on'  t.ituentr;  Pt 

» 

nd 

h- 

We  may  al-o 

;toto  t'- '  na''’-'! 

■  tt^n 

law  in  log  ^it’ii:1  L- 

f  cm 

b i  tehi-'g  log:  'J  t'v.i  of  both  sid>~;  thuc. 


n  K 


r  •  * 

[h  ^  ^ 

-  [x  ’ 


t  I 

+  X  !5o  +  ••• 


n  N-  +  X, 


i  T1  Mr,  +  •  •  • 


(9) 


It  Is  Important  to  note  that  r.  o!  eniool— reoct-ion  equation, 
ri  one  tod  by  (7),ir  simply  a  rc-lntlon  among  the  column  vectors 
Pt  cf  input  elenentr  c.,j  associated  with  each  species  Pj. 

Thus  if  we  let 


“  (aIJ*  a2j *  •••»  ar,j^  * 
wo  have,  corresponding  to  (7), 

X1P1+^Pn+....  X,’  Pj*  +  y  P?*  +  ...  .  (10) 

Tills  is  easily  ween  fron  an  example.  Suppose  H*  and  OlT*  are 
considered  input  "elements”  forming  species  H+,  OH””,  and  H^C. 
We  form  the  tableau 


Input 

S'1  e  front 

Species 

tf+ 

OH 

-J8 

0 

H+ 

i 

0 

1 

X 

Oli 

C 

3 

-j 

where  the  entry  ,  In  row  !  -tnd  column  *  is  the  qumtity  of 
the  input  clement  i  found  in  the  sp  ci.-c  Jf.  It  is  clear  that 
the  chemical  reaction 

H+  +  OH'  HpC 
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is  simply  a  statement  about  the  corresponding  column  vectors 
in  the  tableau,  namely 


We  shall  now  show  that  the  meo3-action  law  holds  fcr  the 
values  of  Xj  that  minimize 

n  _ 

F(X)  -  X  x.F,  ,  (11) 

J-l  J  J 

where  in  particular  Fj  is  given  by 

P.  -  RT  (c.  ♦  In  ^  )  .  (12a) 

J  -lx 

As  noted  earlier,  this  expression  for  Fj  is  essentially  an 
approximation  for  dilute  solutions  and  ideal  gases  of  the  more 
general  expression  (12b).  It  will  be  assumed  unless  otherwise 
stated  that 

Fj  3  RT  ( c j  +  In  a^)  ,  (12b) 

where  each  a^#  the  ’activity,"  is  some  function  of  the  con¬ 
centration  Xj/  x;  that  is  to  say,  a^  (and  hence  F y  the  free 
energy  per  mole  of  the  J-th  constituent)  remains  unchanged  if 
all  quantities  Xj  are  Increased  proportionally.  We  now 
establish  the  following  purely  mathematical  theorem: 


Fj  3  Fj(Xj/  x)  is  a  function 
J=l,  2,  ...,n.  (13) 


4 


’ 
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Because  of  (13)#  F^  is  commonly  called  the  partial  molar 
free  energy  of  the  J— th  CQnrtl tuent .  To  prove  thi3  relation, 
note  first  that  the  free  energy  P(X)  is  a  homogeneous  function 
of  degree  one  in  the  variables  (x^,  x?,  x^),  and  for  any 

homogeneous  functions  of  degree  one  it  is  generally  true  that 


p(x)  =  z  X.  |£-  .  (14) 

>1  3 

This  relationship  (14)  is  known  as  Euler’s  theorem  for  homo¬ 
geneous  forms  of  the  first  degree.  In  words,  it  says  that 
total  free  energy  can  be  determined  a£  if  dF/dXj  is  the 
contribution  of  the  J-th  constituent  per  mole. 

The  foregoing  theorem  is  easy  to  prove;  we  shall  establish 
it,  as  follows:  From  (12b),  the  assumed  homogeneity  of  degree 
xero  for  for  all  t,  we  have 

Xg#  •••#  Xyj)  =  tXg,  •••,  tX^)  • 

Hence,  taking  the  partial  derivative  of  both  sides  with  re¬ 
spect  to  t,  we  obtain 


•  •  • , 


x 


n 


n  ^p(tx^ ,  tXA,  .  • 

)  9  z  - ± ~ - 

1  d<tXj) 


Setting  t  ■  1  yields  (l4). 

To  establish  (13),  we  now  assume  that  F(X)  is  given  by 
(11),  where  P.  ■  F,(x,/  x)  ,  and  take  the  partial  derivative  of  I-'('C) 
with  respect  to  x^.  This  yield; 


i 


* 
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n 


<$F, 


+  w  • 


I  ^  ,  •  •  •  t  !>  I 


(1  ) 


where  we  have  replaced  d?y/dxj  by  d?k/£x  since  *■  1. 

W?  now  show  that  the  summation  t  ;m  vanishes.  Multiplying 
both  sides  by  ;<j  and  summing  on  ,1  yields,  by  relations  ( 1 :* ) 

and  (11), 


0  ■  (x,  +  +  . . .4  x 


"(Jr-'  lr)  • 


Since 


n 

x  >  0,  and  S  x  >  0  , 

•J  ~  ' 

o  * 

the  second  factor  of  the  product  must  vanish  and  (1;)  follows 
from  (15). 

Having  established  that  mdp/fay  let  us  turn  to  the 
main  problem  of  minimizing  the  free— energy  function  F(X) 
subject  to  mas3-balance  constraints, 
n 

^«(X)  *  “  ®n  Xi  -  b,  *0,  i  *  •*■ ,  -#  . .  • ,  m.  .  (1°) 

1  J-l  ’J  1 

Vie  now  review  the  Lagrange  method  for  rini  ilzing  a  general 
function  P(X)  subject  to  r  general  constraints,  (x)  *  0, 
and  later  specialize  the  functions  <J)4  .  The  first  step  is  tc 
assign  unknown  multipliers  ir,  to  the  functions  ♦i  (X)  and  se»v 
the  unconstrained  minimum  of  the  function 


G(X)  -  P(X)  -  »j  (X)  -  li.(X) 


•  •  • 


-  r  %  (X)  .  (17) 
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\ 


Q 

ir  for  some  cnoice  of  m,  ,  the  i  oi..t  C  =  X  whe'*e  the  minimum 

■4. 

of  J(X)  occurs  happens  to  satisfy  the  constraints  ^(X°)  =  0, 
then  this  unconstrained  minimum  point  1°  for  3(X)  is  the  con¬ 
strained  minimum  point  for  ?(X).  To  see  tr.is,  note  that 

ir.in  0(X)  -  G(X°)  •  F(X°)  , 

and  that  for  all  other  X  that  satisfy  |^(X)  =  0  we  have 
?(X)  ^  ?(X°)  because 

F(X)  *  G(X)  £  min  l(X)  =•  ?(X°)  . 

When  F(X)  possesses  partial  derivatives,  and  assumes  its  mini¬ 
mum  at  some  interior  point  of  the  domain  of  definition,  it  can 
be  shown  that  such  multipliers  always  exist.  We  now  specialize 
^(X)  to  be  linear  as  given  by  (l6).  In  this  case,  we  may 
rewrite  G(X)  by  collecting  all  terrn3  with  common  x^,  obtaining 

3(X)  =  F(X)  -  X(ttPj)  *  mb  ,  (lb) 

where 


ir  -  (ir1#  ir2,  . . . ,  ir^) 

is  a  raw  vector, 

^aiJ»  a2J'  amJ^ 

is  the  column  vector  of  coefficients  associated  fith  species; 
J  defined  earlier,  and 

b  -  (fc1,  D;;>  *  *  *'  °n) 

is  the  column  vector  oi  constants.  The  coefficient  of  is 
the  constant  and  dej  anus  on  the  choice  of  m.  At  an  un¬ 
constrained  minimum  jf  'J (  X ) ,  X  satisfies 


for  any  choice  of  v. 


Also  f^r  an-' 


a 


chemical 


'  iven. 


reaction 


(7)  we  may  multiply  the  corresponding  vector  relation  (10) 
by  ir,  obtaining 


X^irP^)  +  X0 ( )  +  ...  ■  A  (irP.,  )  +  A^ '(^o)  +  ••••  (19) 


which  holds  for  any  choice  of  v,  Tims  substituting 
dF/dXj  ■  tPj,  we  have 


$E_ 

6t, 


+  A 


dF 

6c, 


•  •  •  * 


(20a) 


which  holds  for  the  X  that  minimizes  G(X)  for  any  choice  of 

tt.  Since  P  has  a  special  form  (11)  that  satisfies  dP/dx*  ■ 

J  J 

we  also  have,  for  any  such  X, 

X1  *1  +  XJ  *2  +  "•  “  Xl'*i  +  W  +  •"  •  .  {50b) 

If  we  further  set 


Fj  -  RT  (c  T  +  In  cij )  , 

then  upon  substitution,  rearranging  of  terms,  and  dropping  of 
tiie  common  factor  RT,  we  get  the  relationship 

•  t  it 

(A-,c^  +  XpC2  +  •••)  —  (A.,C|  ^  ■t  •••) 

(1C 

i  ii  i 

■  (Aj  In  In  +  ...)  -  (A^  In  In  a,,  +  .,.) 

for  all  X  tnat  minimize  Q ( X )  for  any  choice  of  w.  In  par¬ 
ticular,  for  that  tt  that  yields  X  -  X°  satisfying  the  moss 
balances  4>^(x)  ■  C  and  thereby  minimizing  the  fr^e  energy 
F(X),  this  sane  relationship  must  hold.  But  this  holds  tven 
if  we  start  out* with  a  different  set  of  b^  values  since  the 
expression  is  independent  of  b, ,  If  w;  'ine  In  K  to  be 


•  •  *  )  »  (2-) 


II  II 

In  K  ~  ^  X ^0 ^  •  •  •  ]  “  |X«  ^  X/>C/^ 

X  -i  d  c.  IX  d  *_ 

then  (20c)  establishes  the  mass-action  law  In  logarithmic 
form  and  In  terns  of  activities.  If  the  activities  are  given 
by  ■  x^/  x,  then  this  is  the  sane  as  the  mass-action  law 
stated  earlier  in  equation  (9)«  If 

aj  ■  V  * » 

where  is  some  constant,  this  also  yields  the  mass-action 
law  stated  earlier  with  an  adjusted  value  for  In  Ki.  In  this 
manner,  the  values  from  the  free-energy  function  can  be 
used  to  define  the  equilibrium  constants.  The  converse  is 
shown  below. 

DETERMINING  Cj  VALUES  FROM  EQUILIBRIUM  CONSTANTS 

As  was  noted  in  the  text,  the  equilibrium  constants  for 
many  common  reactions  are  tabulated  in  physical  chemical  tables. 

We  wish  to  discuss  how  to  adjust  the  Cj  values  so  that 
the  equilibrium  constants  are  directly  applicable  to  the  model;- 
to  be  precise,  we  note  that  for  those  x  that  satisfy  the  mass 
balances,  the  expression  for  F(X)  may  be  replaced  by 

P(X)  -  "  x,  ( c .  +  m  j)  -  "  ::,b,  ,  (22) 

J-l  3  1  i-1  1  1 

where  the  k^  are  constants  and  the  bear  a  direct  relation¬ 
ship  to  the  equilibrium  constants.  The  numbers  k^  wer*-* 
selected  so  that  when  multiplied  by  the  i— th  mass-balance 
mtion  and  the  m  equations  summed  and  subtractec  from  F(X). 


it  of  the  new  coefficients  of  the  x^'s,  denoted 
vanish.  In  place  of  the  original  function  P(X) 


c y  would 
,  the  function 


o(x)  -  "  x,(3,  +  in  ii)  (?.') 

J-l  1  3  * 

was  minimised  subject  to  the  mass-balance  constraints.  Since 

the  two  functions  differ  by  a  constant, 

m 

2  k.b.  , 

1-1  1  1 

the  value  of  X  -X0  that  minimizes  P(x)  also  minimises  G(x). 

Any  m  of  the  cT^  can  be  made  to  vanish  by  suitably 
choosing  the  values,  provided  that  the  square  array  of  co¬ 
efficients  | ^ |  associated  with  these  J  Is  nonsingular,  i.e., 
provided  that  their  determinant  /  0. 

Por  Model  I,  the  Oj  values  associated  with  species  H+ 
and  0H"~  in  Compartment  II  were  made  among  others  to  vanish. 


CH+  *  C0H~  "  0  • 

Now  the  equilibrium  constant  for  water  does  not  depend  on 
whether  we  use  c1  or  c .  values  for  Its  definition  in  (21), 

u  •) 

Its  value  being  the  same  since  it  must  yield  the  same 
right-hand  member  of  (20c).  Hence,  for  Compartment  II  ipecles. 


ln  KH20  *  cHo0 


CH+  “  C0H  "  CH?0  ‘ 


In  other  words,  the  value  of  c\  0  «  cu  „ 
obtained  from  the  equilibrium  constant 
C,  the  entire  set  of  n  —  m  remaining 
in  this  manner  from  equilibrium  conot^*) 


•an  be  directly 

or  water.  In  Appendix 

values  is  evaluated 


No  attempt  was  made  In  the  construction  of  Model  I  to 
determine  the  absolute  c.  values  as  defined  in  terms  of  Ff 
value 3 .  All  Cj  values  shoirn  are  the  relative  values  as 
obtained  from  equilibrium  constants . 


i 
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Appendix  B 

CONSTRUCTION  OF  THE  MODEL  —  INPUTS 

The  theoretical  background  given  in  Appendix  A  makes 
possible  a  more  detailed  description  of  the  construction  and 
operation  of  the  model  shown  in  Table  1.  The  first  block  on 
the  left,  labeled  "Inputs/’  contains  the  "elements"  from  which 
the  molecular  species  in  the  three  compartments  —  "Air  Out," 
"Arterial  Plasma,"  and  "Red  Cells"  —  are  constructed.  These 
Input  "elements"  have  the  values  in  moles,  b1#  b11#  that 

are  present  per  unit  of  time  in  mixed-venous  blood,*  plus  the 
values  in  moles  of  fresh-air  gases,  a^,  ...,  ac,  entering  the 
lung  sacs  per  unit  of  time.  For  the  validation  of  the  model, 
one  liter  of  average  mixed-venous  blood  was  used  as  an  input, 
plus  ~  0.031  moles  of  degraded  "fresh"  air. 

Inputs  bj,  b2,  b^,  a^,  a2,  a y  are  obvious.  The  input 
"elements"  in  rows  4  and  5,  H+  and  0H~,  provide  a  convenient 
way  to  form  water  in  the  several  compartments  and  serve  as  a 


A  troublesome  data  problem  arose  in  attempting  to 
validate  the  model.  Even  among  normal  individuals,  the  varia¬ 
bility  of  blood  composition  is  great,  and  in  addition  it  is 
necessary  for  validation  to  have  laboratory  data  for  venous 
and  arterial  blood  on  a  completely  comparable  basis. 
Fortunately  for  our  purpose,  the  U.S.  Air  Force  has  prepared 
compilations  of  biological  data  in  this  form.  The  data  for 
Inputs  and  for  comparison  with  outputs  were  obtained  from 
Refs.  9  and  13. 

•  • 

As  explained  elsewhere,  the  air  inputs  may  also  be 
obtained  by  using  a  large  amount  of  air  of  lung-sac  composition 
or  by  programming  the  computer  to  search  for  the  proper  amount 
of  "fresh*  air  to  meet  blood  or  lung  criteria. 
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mass  source  and  sink  for  the  reactions  Involving  H+  and  OH”, 
for  example, 

C02  +  OH”  HCO”  . 

It  will  be  noted  that  HCO”  Is  not  one  of  the  Input  "elements," 
although  It  Is  an  Important  constituent  of  blood.  The  model 
creates  this  species,  however,  through  the  reaction  given 
above.  To  provide  for  this  HCO^  formation  in  terms  of  mass 
and  charge,  the  input  values  of  H+  and  OH”  are  set  different 
by  the  amount  of  HCO”  that  is  present  in  mixed— venous  blood. 
The  C02  input  b2,  which  Is  also  involved  in  production  of 
HCO”,  is  the  total  C02,  including  that  present  in  combined 
fona.  It  i3  not  directly  involved,  of  course,  in  the  charge 
balance.  The  mass-difference  effects  in  terms  of  H*  and  OH” 
caused  by  the  slight  basicity  of  the  blood  is  so  small, 
approximately  1  x  10”10  moles,  that  it  does  not  require 
separate  mas3  accounting  at  this  stage. 

"Elements"  Na+  and  K+,  in  rows  7  and  8,  are  introduced 
at  present  essentially  as  fixed  values  in  plasma  and  red 
cells,  respectively,  inasmuch  as  entries  for  these  species 
occur  only  and  separately  in  these  two  compartments,  and  are 
not  allowed  to  migrate  in  the  model.  In  real  life,  of  course, 
there  is  some  K+  in  plasma  and  Na+  in  red  cells.  There  are 
also  other  cations  in  both  compartments,  as  established  by 
the  selective  screening  action  of  the  red-cell  membrane.  If 
and  when  the  data  may  be  available,  this  lon-ecreenlng 


•1 


tr 


l 


mechanism  of  the  membrane  can  be  incorporated  in  the  model  as 
a  series  of  chemica)  equations  with  their  accompanying  thermo¬ 
dynamic  functions.  At  the  present  stage,  however,  all  cations 
(except  H+)  in  plasma  and  red  cells  are  stated  in  terms  of 
equivalent  Na+  and  K+,  respectively. 

The  "element"  in  row  9#  designated  by  the  shorthand  HBl“, 

represents  the  total  hemoglobin  present  in  red  cells.  It  is 

# 

introduced  at  this  stage  as  moles  of  one-fourth  hemoglobin 
molecule,  molecular  weight  approximately  16,300  (one  02, 
one  Fe  per  unit).^  The  present  model  does  not  allow  for 
the  hydrogen  ionization  of  hemoglobin.  The  aggregating 
techniques  suggested  by  Alberty^)  provide  a  method  for  in¬ 
corporating  the  multistaged  ionization  of  hemoglobin  Into  the 
model  with  the  addition  of  a  modest  number  of  columns. 
Unfortunately,  the  data  provided  for  the  several  ionization 
constants  are  for  horse  hemoglobin  at  23°C.  If  similar  data 
for  the  ionization  of  human  hemoglobin  become?  available,  it 
can  be  incorporated  with  possibly  some  Improvement  in  the 
representation  of  "acidity."  In  Mode]  II  (see  Appendix  D) , 
we  have  followed  an  approach  suggest  y  Pauling^1-^  with 
regard  to  the  ionization  of  the  hemoglobin  molecule. 

The  "elements"  in  rows  10  and  11,  protein  in  plasma  and 
In  red  cells  (other  than  hemoglobin),  are  Inserted  in  the 
model  at  this  stage  in  terms  of  moles  of  average  molecular 
weight  (equivalent  moles)  to  represent  the  wide  range  of 
molecular  species  of  protein  that  actually  occur  in  the  b^ood. 


At  a  later  stage,  and  when  the  data  permit,  these  protein 
species  may  be  Introduced  separately  and  their  Ionization  may 
be  represented  in  the  same  fashion  a 3  suggested  for  the  case 
of  hemoglobin. 

As  mentioned,  row  12  contains  the  charge  equation,  which 
represents  the  conservation  of  charge  and  accounts  for  the 
Olbbs-Donnan  equilibrium  across  the  red-cell  membrane.  At 
present  we  have  taken  b12  «=  0,  which  dictates  electrical 
neutrality  among  the  several  compartments.  Any  charge  dif¬ 
ference  that  may  occur,  as  between  red  cells  and  plasma,  can 
be  accounted  for  and  its  influence  on  relative  composition 
represented  by  the  Introduction  of  the  appropriate  molal 
charge-difference  value  of  b12.  Charge  differences  of  the 
order  of  -16.8  millivolts  between  red  cells  and  plasma  are 
reported.  Charge  differences  of  this  magnitude  could  have 

a  measurable  and  possibly  physiologically  significant  effect 
on  the  composition  of  plasma  and  red  cells.  Exploration  of 
this  phenomenon  by  means  of  the  model  remains  for  the  future. 
With  b12  *  0,  the  chloride  shift*  occurs,  so  that 

plasma  _  feared  cells 

frco>] plasma  ^HC05red  cells 


#It  can  be  shown  that  thl3  same  relation  holds  if  we 
allow  the  net  charge  b12  in  Compartment  II  to  be  variable  and 

assume  that  the  energy  required  to  maintain  this  charge  can 
be  expressed  by  adding  a  new  tenn  to  the  free-energy  function. 
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Appendix  C 

DETERMINATION  OP  FREE-ENEROY  VALUES 


So  far,  we  have  described  the  operation  of  the  model  in 
terms  of  the  several  compartments,  the  thermodynamics  involved, 
the  inputs,  and  the  mass  and  charge  conservation  constraints. 

We  shall  now  discuss  the  numerous  chemical  reactions  involved, 
how  they  are  represented  in  the  model,  and  how  the  mass-action 
or  partlal-mola!  free-energy  constraint  operates  (with  the 
others)  to  yield  an  equilibrium  solution.  The  mass-action 
equations  describe  the  formation  of  the  molecular  species  of 
columns  .  .  .  ,  in  three  compartments,  I:  Air  Out.. 
II:  Arterial  Plasma,  and  III:  Red  Cells,  and  relate  similar 
species  occurring  in  more  than  one  compartment. 

For  example,  the  entry  1  in  column  12,  row  4,  Compartment 
II,  the  entry  1  In  the  same  column,  row  5*  and  the  entry  c^ 
in  the  bottom  row  together  express  the  chemical  equation, 

H+  +  ON"  ^  H20  , 


In  the  mass-action  form, 


In 


tH2°hl 


°12  * 


where 

"N 

Ci2  “  ln  *W' 

and  is  the  ionization  constant  for  water  (on  the  mole- 
fraction  scale)  at  >7°C.  It  was  shown  previously  that 


< 


instead  of  the  repara te  relations 


plasma 

red  cells  * 

and 

tHC0P  plasma 

[HC05  red  cells  * 

which  would  hold  if  the  charge 

equation  were  removed. 

The  input  values  of  v1#  . 

*  *  *  v12 '  al'  *  *  #  *  aS» 

and  t>lt 

1 


,  b12,  for  which  the  equilibrium-solution  values 
30 


were  obtained,  are  listed  in  Table  2  in  the 


text. 


P-1811 

71 


constants  like  c12  derived  from  equilibrium  constants  may  be 
used  In  the  model  with  the  3ame  results  as  if  the  individual 
F°/(RT)  values  for  each  species  were  obtained  from  tables  and 
used. 

# 

In  a  manner  similar  to  that  described  for  the  water  re¬ 
actions  in  Appendix  A,  the  other  reactions  involved  are 
entered  in  the  model.  Por  the  same  species  occurring  in  both 
Compartments  II  and  III,  the  assumption  ir  made  that  their 
original  energy  coefficients  (not  necessarily  c,  values) 
are  the  same  in  both  compartments.  Thi3  results,  for  example, 
in 

[°2]  ii  *  [°2]  hi  • 

The  other  reactions  and  their  Cj  values  are  given  in  Table  4. 
In  this  table,  the  species  are  shown  in  square  brackets,  [  ], 
indicating  concentration  in  mole  fraction.  The  c^  values  were 
calculated  from  the  K  values  shown  in  Table  5,  all  converted 
to  the  mole-fraction  scale. 

It  Is  Important,  in  applying  these  data  to  the  model, 
that  the  K  values  from  the  literature  all  be  converted  to  the 
same  scale  (i.e.,  molar,  molai,  or  mole  fraction).  Most  K 
values  are  reported  on  the  molal  scale,  some  on  the  molar; 
and  In  some  cases,  unfortunately,  the  scale  is  not  indicated. 
All  K  values  were  converted  to  the  mole— fraction  scale  for 
U3e  In  the  model.  It  is  also  Important  that  the  K  values  are 
those  for  the  temperature  of  interest  (^10.18  Kelvin  in  the 
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Table  k 

REACTIONS  AND  FREE  ENERGY  VALUES a  POR  MODEL  I,  Z7°C 


Column*1 

Mas3  Action 

=  ,  -  m  k 

(i) 

In  IX  /  ^2^1 

or  III 

-10.  So 

(2) 

In  [CO,]  n  /  [COjj 
or  III 

i 

-7 .  *.'9 

(?)  | 

In  [Np3  j j  /  [N^] i 

or  III 

-r 

| 

(12) (2*0 

In  /  [H20]  XI 

or  III  or  III  or  III 

-50.2? 

1 

(*)  : 

In  [H+] n  [OH-]n  /  [H,0] t 
or  III  or  III 

j 

1 

* 

-  c,2  +  In  [«20]n  /  [HgOlj 
or  III 

-56.44 

(1?) (25) 

1 

| 

in  [C02]  „  [ort-]  I:  /  [hco-]jt 

or  III  or  III  or  III 

-21.20 

(29) 

In  [Ojj] j jj  /  [HBlOglm 

1 

•  +  In  [Og]  j  [HB1  J  Hi  /  [HBlOjx]  jjj  | 

-16.25 

(1")(26) 

In  [HgOj  j  j  / 

or  III  or  III  or  III 

0 

05)  (27) 

_ 

In  II  /  [coj] ii  ji 

or  III  or  III  or  III  or  III 

i  ..  ■■■  —  ■  —  ■  —  ■■  -  — . 

g.2*: 

All  Cj  values  shown  are  the  adjusted 

to  In  Appendix  A. 


values  referred 


bThe  valued  of  c.  ■  0  for  columns  (  ),  {'  ),  ("),  (8), 
(?),  (10),  (11),  (I?),  (17),  (IS),  (19),  (20),  (21),  {?■>), 
(25),  (28),  (JO). 


< 
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Table  5 

EQUILIBRIUM  CONSTANTS* 


Column 

Equilibrium  Constant 

0) 

K  -  0.0232b 

(2) 

K  -  0.5672° 

(5) 

K  -  0.0127* 

(12) (24) 

l°e10  K  -  -lj.5bd 

(“t) 

l°fi>10  ^  *  ^°®10  *w  ^  *°®10  II  /  OV  • 

or  III 

log10  /  [HgOlj  -  1.2U96e 

or  III 

(U)  (25) 

log10  K  -  log10  1C*  -  log10  IC^  co  . 

log10  K.  ■  c.09f 

HgCO^ 

(29) 

log10  K  -  -7.0486® 

( 14 ) (26) 

log10  X  -  700h 

(15) (27) 

lo®10  K  ■  log10  ^C02^II  ^0ir^II  /  ^HC05h^K2h  co  ^ 

2  3 

log,n  K9  -  -10. 251 

10  %co5 

aAll  K  values  on  the  molar  or  molal  scale  were  converted  to 
the  mole  fraction  scale  for  calculating  c's. 


brtef .  3,  p.  54. 
cRef .  16,  p.  1092. 
dRef.  17,  p.  152. 
eRef .  16,  p.  1465. 
fRef.  9,  p.  27 2. 
®Ref.  3,  p.  64. 
hHef.  3,  P.  80. 
1hef.  16,  p.  1199. 
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model  30  far),  or  are  converted  to  the  value  for  K.^  1g  If 
the  necessary  thermodynamic  data  permit.  At  a  later  stage, 
the  appropriate  thermodynamic  functions  can  be  introduced  into 
the  model,  relating  free  energy  to  temperature  and  pressure, 
and  perhaps  to  electrical  and  gravitational  fields  and  to 
surface  effects. 


Mention  was  made  in  Sec.  VII  of  certain  deficiencies  in  our 
i^odei  I  of  the  respiratory  function.  those  relating,  to  the 

chemistry  of  the  tlood,  the  most  important  are  the  representa¬ 
tion  of  hemoglobin,  including  the  carbamino  reaction,  an^  the 
mul  tistaged  ionization  of  protein  and  protein— 1  Ike  substances 
in  the  plasna  and  red  cells.'  An  expanded  vodel  II,  shown  in 
Table  6,  was  constructed  to  demonstrate  one  of  several  ways  in 

• 

which  the  representation  of  these  phenomena  might  be  improved. 

In  this  model,  two— staged  ionization  of  proteins  is  ac¬ 
counted  for  In  both  plasma  and  red  cells.  The  representation 
of  this  type  of  phenomenon  could  be  expanded  to  any  degree  of 
detail  desired.  In  addition  to  proteins,  it  might  include  a 
wide  range  of  individual  molecular  species  such  as  amino  acims; 
and  it  might  include  other  reactions  than  ionization  —  e.g., 
esterification. 

In  Model  II,  the  more  detailed  structuring  of  nemoglooin 
and  its  reactions  is  represente i  by  a  separate  compartment  IV, 
as  if  a  separate  phase  were  involved.  There  ire  distinct  ad¬ 
vantages  to  this  form  of  representation  for  such  complex 
structures  and  reactions  whenever  it  can  oe  assumed  that  re¬ 
actions  with  different  parts  of  the  molecule  operating 


\ 


t 


X* 


r, 

j 

.■* 

■I  r 


> 

4. 


,  '* 
i 

T> 


A  numoer  jV  hyj  obliges  xist  In  the  I  i  teraturc-;  /.e  have 
cnosen  those  due  to  Adair  (•)»  Paulin,  { 1  ),  iaa  -ovghton 
among  others,  for  1  iustrut 1 v-  |  arm  ;ses  o n'v.  In  cases  wriepe 
there  are  competing  hyj  ot.heses ,  a  ■  terr.at  Ive  ~  o  .jl  c  >ul  i  c  ^ot 
u;  and  used  t  see  %-nic.n  test  j  redl.cts  co^r  ex  ser1  ed 

rh  no  err  . 


>DEL  II  OF  THE  EXTERNAL  RESPIRATORY  oYSTEM 
OUTPUT  COMPOSITION  BY  MOLECULAR  SPECIES 
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4 


05  O 
1  : 

• 

•  • 

•  • 


refers  to  an  a  Id,  nonoxylablle  group  of  hemoglobin, 
refer  to  an  oxygen— 'onnectlng  group  cf  a  one— fourth  hemoglobin  koi«:ule. 
refers  to  gn  acid,  oxylablle  group  cf  hemoglobin. 
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independently.41  Formal  compartmentalization  is  not 
essential  —  by  a  somewhat  different  arrangement,  hemoglobin 
and  its  reactions  could  be  incorporated  in  Compartment  III, 

Red  Cells.  Columns  24—27  and  39—44  were  designated  aa  species 
of  hemoglobin,  and  as  portions  of  its  molecule,  to  correspond 
with  Adair's  hypothesis^)  of  four  levels  of  oxidation  of  the 
hemoglobin  molecule  and  with  Linus  Pauling'3  ;qaare-heme- 
3tructure  hypothesis  and  techniques  for  accounting  for  the 
variable  affinity  of  oxygen  with  changes  in  oxygen  concentre-'' — 
tion  and  with  changes  in  pH. 

Column*  }8,  and  }4  account  for  the  CO^— carbamino  link¬ 
ing  with  hemoglobin  and  its  relation  with  oxygen  and  hydrogen- 
ion  concentration.  The  mechanisms,  puggested  by  Houghton 
(pp.  86-88  of  Ref.  4)  have  been  followed  in  this  representation. 

These  more  detailed  reactions  of  hemoglobin  as  repre¬ 
sented  in  compartment  IV  may  be  described  briefly  as  follows: 
Ihe  six  possible  molecular  arrangements  of  hemoglobin  accord¬ 
ing  to  Pauling  are  here  designated  as  "B",  0^"B2"(a), 

0^MB2"(b),  and  (columns  29—44).^  For  those 

molecular  arrangements  (columns  4l,  4j>,  44)  where  the  square 
oxygen— linked  heme  groups  Interact  nonaymetrically,  the  free 

#Without  cu^h  Zi  devi  e.  It  would  be  neofcssary  to  set  up 
separate  columns  to  represent  perhaps  a  million  distinct 
variants  of  the  hemoglobin  molecule  corresponding  to  the 
different  ways  its  more  than  ICO  groups  and  parte  can  ionise, 
or  bond  with  C0c  or  0^.  ’ 


P-1811 

78 


energy  is  decreased  algebraically  by  a  factor,  RT  In  a,  for 
each  unbalanced  aide,  where  a  is  an  empirical  value  determined 
from  experimental  evidence.  Inasmuch  aa  AF  ■  -  RT  In  K  , 

i 

where  K  i3  the  free-energy  change  accompanying  the  addition 
of  oxygen  to  heme,  the  "c"  values  for  the  model  may  be  ob¬ 
tained  from  the  logarithms  of  the  separate  and  appropriate  a 
factors  for  each  of  the  molecular  arrangements. 

an  analogous  manner,  Pauling  suggests  that  the  effects 
of  .ges  in  hydrogen-ion  concentration  on  hemoglobin  may  be 

i 

effected  by  an  empirical  correction  to  the  K  value.  In  this 
case,  empirical  factors  RT  In  0  and  a  30-called  acid— strength 
constant,  RT  In  A,  are  used  to  account  for  the  difference  in 
hydrogen-ion  dissociation  of  oxygenated  and  unoxygenated 
hemoglobin.  These  correction  factors  are  determined  from 
experimental  evidence.  The  precise  mechanism  for  doing  this 
is  discussed  below. 

In  terms  of  Model  II,  columns  24  and  22#  "A~"  and  H'A" , 
represent  the  groups  of  the  hemoglobin  molecule  that  can 
lonice  and  that  are  not  affected  by  oxygen  bonding.  There 
are  n ^  of  these  hydrogen-bonding  groups.  The  energy  constant 
per  group  is  assumed  to  be  Pauling's  In  A.  Columns  2o  and 
27 #  "C“"  and  H"C",  represent  an  additional  eight  ionizing 
groups  of  the  hemoglobin  molecule  near  to  and  affected  by 
oxygen  bonding  ( "oxylablle") .  Two  acid  groups  interacting 
with  each  heme  were  postulated  by  Pauling  to  obtain  satis¬ 
factory  agreement  between  the  model  and  experimental  results. 
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Each  "oxylabile"  group  is  assumed  also  to  have  the  energy 
constant  In  A  when  no  oxygen  is  attached  on  the  heme  and 
In  A  -  In  3  when  oxygen  is  present.  Following  Houghton,  we 
allow  one  COp  to  attach  to  one  of  two  oxylabile  group*  If  .10 
Op  is  present. 

•  1 

From  the  empirically  determined  values  for  K  ,  a.  A,  and 
0  from  Pauling,  It  is  therefore  possible  to  calculate  the 
values  -  c^  for  use  In  Compartment  IV  of  Model  II.  The 
mass  balance  of  all  the  hemoglobin-associated  species  of 
Compartment  IV  is  related  to  total  hemoglobin  in  Compartment 
III,  Red  Cells,  through  column  33  of  Red  Cells  and  through 
the  Op,  COp,  and  H*  coefficient  entries  of  the  hemoglobin- 
associated  Input  "Elements." 

Model  II  has  not  been  validated;  until  this  is  done  it 

1 

cannot  be  certain  that  K  ,  a.  A,  and  3  values  have  been 
correctly  interpreted  in  the  model.  For  example,  the  energy 
constants  of  Model  II  are  on  a  per-group  basis,  whereas 
Pauling's  A  and  3  values  may  apply  to  the  computation  of  a 
pair  of  groups. 


I 


i 
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Appendix  E 

SIMULATING  THE  RESPIRATORY  FUNCTION  ON  AN  ANALOG  COMPUTER 
THE  USES  OF  THE  COMPUTERS 

Computer  solutions  of  the  lung  model  have  been  undertaken 
with  two  goal  a  In  mind.  The  first  Is  the  obvious  necessity 
to  make  sensitive  checks  on  the  mathematical  model.  A 
physiological  model  13  necessarily  founded  on  many  parameters 
and  a  number  of  hypotheses,  both  explicit  and  implicit.  The 
ability  to  make  fast,  deliberate  tests  of  these  underlying 
assumptions  has  aided  materially  in  the  validation  of  this 
first  model. 

The  second  goal  is  to  obtain,  in  physical  hardware,  a 
simulation  of  the  time-dependent  system.  Such  a  simulation 
should  present  the  complex  interplay  of  component  chemical 
reactions;  but  it  should  do  more.  It  should  be  amenable  to 
the  analysis  of  individual  reactions  as  well  as  the  whole; 
it  should,  of  course,  compute  the  equilibrium  conditions;  and 
it  should  have  sufficient  capacity  and  spe-d  to  enable  an 
analysis  of  the  transient  change  from  one  equilibrium  state 
to  another.  Such  a  change  in  equilibrium  conditions  will 
occur,  for  example,  when  the  lungs  are  suddenly  exposed  to  pure 
oxygen  instead  of  the  normal  air  mixture. 

Almost  any  physiological  model  that  approaches  reality 
will  be  sufficiently  complex  to  obfuscate  one's  mental  image 
of  the  3y3tem.  But  the  analog— computer  system  has  the 
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advantage  that  it  may  be  built  up  a  piece  at  a  time  to  any 
desired  complexity.  Hence,  the  goals  mentioned  above  seem 
perfectly  feasible,  although  certain  limitations,  such  a3  com¬ 
puter  size,  are  apparent.  Despite  the  fact  that  the  two  goals 
overlap,  it  was  decided  to  attain  them  on  different  computers. 
The  digital  computer,  with  high  ac  uracy  but  relatively  long 
solution  time,  is  used  to  compute  the  equilibrium  conditions 
of  the  system  and  hence  to  test  the  underlying  hypothesis  of 
the  model.  TTie  analog  computer,  with  lower  accuracy  but  much 
shorter  solution  time,  can  be  used  to  search  through  a  range 
of  parameters  or  ‘an  eventually  be  used,  we  believe,  to  simu¬ 
late  the  dynamic  aspects  of  the  model.  Por  the  present,  it 
is  a  valuable  tool  for  tracing  the  effects  of  continuously 
varying  inputs  to  the  system.  Por  example,  the  analog  can  be 
U3ed  in  complex  experimental  situations  for  estimating  equi¬ 
librium  constants  by  regarding  these  constants  as  parameter.1 
to  be  varied  continuously  cn  the  machine  until  a  good  fit  is 
found  between  observations  and  predictions. 

CONTINUOUS  OPTIMIZATION  AND  THE  ANALOG  COMPUTER 

The  problem  of  continuous  optimization,  as  compared  to 
the  static,  one-time,  optimization  of  Appendix  D,  is  to  com¬ 
pute  and  display  the  transient  phenomena  of  a  complex  system 
as  it  changes  under  the  influence  of  a  time-dependent 
parameter.  The  central  postulate  is  Hamilton's  principle, 
the  principle  of  least  action  from  classical  mechanics. 
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The  system  Is  defined  as  In  the  text.  The  free-energy 
function  P(x)  Is  to  be  minimized  subject  to  mass-balance  and 
charge  constraints.  Using  the  Lagrangian,  we  write,  for  k  «  I, 


IT,  III, 

0(X,t)  -  P( X)  -  ",  T.$,  -  2  x  (c ,+ln  -  I  r,4>.(x),  (24) 

1-1  1  J-l  J  J  xk  i.i  1  1 


where  ♦t(x) 


n 

Z  a.  .x.  -  b. .  At  equilibrium,  we  require 
J-l  J  1 


■  c.  +  In  -  2  a. .t.  -  0  for  each  J, 
J  \  1-1  1 


(25) 


-  ♦,  (X)  -0  for  each  1, 

*1  J 

where  the  are  the  Lagrangian  multipliers. 

A  unique  solution,  or  equilibrium  condition,  obtains  on 
the  computer  when  a  set  of  mole  numbers  X  and  a  vector  v 
Identically  satisfy  the  equations  (25). 

If  the  physical  system  represented  by  Model  I  Is  well 
behaved.  It  will  exhibit  at  a  solution  certain  stability 
properties,  and  the  system  will  remain  stationary  until  a 
parameter  of  the  system  13  redefined,  e.g.,  the  quantity  of 
a  reactant.  At  such  a  time,  transient  phenomena  are  exhibited 
until  a  new  equilibrium  obtains.  Por  most  parameters  of  the 
system,  small  changes  of  the  parameter  reflect  small  changes 
in  the  solution;  that  is,  the  solution  Is  a  continuous  func¬ 
tion  of  the  parameters  of  varlatio  1.  Thud,  if  a  parameter  is 


a  slowly  changing  function  of  time,  the  solution  moves  In  a 
trajectory  that  continuously  minimizes  F(X).  Thus  (25)  always 
holds,  although  0  *  0(X,ir,t). 

THE  TECHNIQUE  OF  THE  ANALOG  SOLUTION 

Accurate  continuous  simulation  of  a  :hemi  al  complex 
may  be  accomplished  on  an  analog  computer  by  means, 
essentially,  of  the  method  of  steepest  descent.  The  analog 
Is  desirable  because  of  the  characteristically  short  solution 
time  available,  particularly  in  the  high-gain  mode  of  •ompc- 
tation.  The  general  method  is  shown  in  a  previous  paper, 
"Continuous  Programming  Methods  on  the  Analog,  which 

contains  a  bibliography.  Some  modifications  have  ensued, 
however,  which  represent  a  departure  from  previous  pro¬ 
gramming  methods;  some  details  will  accordingly  “be  shown  in 
this  Appendix. 

The  essential  Justification  for  the  method  employed  come 
from  the  studies  of  convergence  questions  in  the  litera¬ 
ture.  (^9) (20)  ^  Kose,  in  a  1956  paper,  considered  the 

convergence  of  a  solution  of  the  mathematical  programming 
problem  as  solved  on  an  analog  computer.  To  use 
his  result,  it  is  necessary  to  postulate  the  existence  and 
uniqueness  of  a  solution,  considerations  whi’h  are  taken  up 
for  more  general  programming  problems  elsewhere  in  the 
literature  of  mathematical  programming.  With  these 
hypotheses,  we  have  the  following  result. 


P-1811 

84 


Theorem.  If  G(X,t)  Is  a  strictly  concave,  dl  ’arentlable 
function  of  X  as  well  as  a  strictly  convex  and  dlfi ^rentlable 
function  of  r  for  >  0,  then  all  solutions  of  the  system  of 
of  differential  equations. 


a 


x 


~a  t 


1 


da 


5 


6,  i(x)  , 

Ti  1 


it  Xj  m  0  and 
otherwise  , 


da 

dXj 


do 


<  0  , 


Iff.  «=  0  and  ^-  >  0  , 

&▼! 


otherwise  , 


a  «  a  positive  time  constant, 
necessarily  converge  to  a  unique  saddle  point. 

It  can  be  shown  that  the  conditions  of  the  theorem  are 
satisfied  by  the  function  0.  Ihe  constant  a  Is  arbitrary  and 
we  wish  to  make  It  very  small  so  that  the  solution  may  be 
obtained  In  a  short  time.  Of  course,  various  practical  con¬ 
siderations  obtain  on  the  analog  computer,  because  the  com¬ 
puter  Is  a  physical  device  with  characteristics  (e.g,, 
frequency  limitations)  of  Its  own,  and  these  characteristics 
supervene  on  the  theoretical  mathematical  statement.  With 
these  considerations  In  mind,  we  may  proceed  to  program  the 
problem. 


t 

i 
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It  Is  required  to  maintain  the  relations  (25)  Identically, 
independent  of  a  new  time  variable  t.  To  this  end,  the  rela¬ 
tions  (25)  are  programmed  on  amplifiers  in  the  high-gain  mode, 
the  outputs  of  which  are  the  vectors  X  and  ir.  In  order  to 
explain  this  mode,  we  start  with  the  basic  operational  unit 
of  the  computer,  the  high— gain  amplifier,  represented  by  the 
symbol 


where  is  the  "gain"  of  the  amplifier  and  is  usually  of  the 
order  of  l()f .  Note  the  sign  change  and  the  addition.  Sin  e 
the  output  variable,  e^,  is  limited,  for  practl-al  ele:troni 
reasons,  to  +  100  units,  the  sum  of  the  inputs  is  bounded  by 
±  10  units;  l.o.,  the  sum  of  the  inputs  must  be  very  near 
zero,  and  we  regard  this  device  .  s  a  rectrl  tion  on  the  In¬ 
puts  . 

In  order  to  control  the  sura  of  the  inputs,  the  amplifier 
is  normally  used  with  a  "negative  feedback"  loop.  The  design 
of  this  loop  defines  the  mathematical  operation  of  the  unit. 
Thus,  Just  as  an  example,  it  _s  possible  to  take  the 
(positive)  square  root  of  an  input,  e^,  by  inserting  a 
"multiplier"  in  the  feedback: 
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A  sign  change  occurs  In  the  multiplier.  A  heuristic  explana¬ 
tion  of  the  operation  may  be  had  by  regarding  the  amplifier 
a  an  integrator.  The  output,  i.e.,  the  Integral,  will  be 
steady  when  the  sum  of  the  inputs,  i.e.,  the  derivative,  ir- 
zero.  Thus,  A ( — x  tx  )  Is  a  differential  operator  for  e0, 
the  whole  system  is  a  function  of  time,  and  as  t  gets  large 
A(-x  +x*" )  goes  to  zero.  This  system  is  stable  since  if  ej^ 

is  given  a  positive  increment,  e^  takes  a  negative  increment 

2 

and  -x  de  rease:  bringing  the  sum  of  the  inputs  again  to 
zero.  The  whole  operation  should  be  regarded,  now,  as 
determining  e^.  If  the  gain  i3  sufficiently  high,  say  10  , 
the  accuracy  is  very  good  since  small  errors  in  cause 

very  significant  errors  in  the  zero  sum  of  the  inputs  and 
in  the  right  direction  to  correct  e^. 

An  Important  thing  to  note  is  that  was  arbitrarily 
defined;  it  has  no  necessary  relation  to  e^.  But  once  it  is 
defined,  =  f{«i),  it  is  necessary  to  construct  the  inverse 
fui  ion  generator,  e^  »  f  (eQ;-  for  the  feedback  loop. 

Also  the  feedback  loop  may  itself  be  implicitly  defined. 

Using  such  self-regulating  loop-,  we  may  in  trument 
the  complex  equilibrium  problem,  and  the  'omputer  will 
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automatically  seek  the  equilibrium  valuer.,  i.e.,  values  for 
which  each  of  the  Interconnected  loops  is  satisfied.  Such  a 
solution  requires  a  few  milliseconds.  Thus,  we  instrument 
the  mass-fcalance  equations, 

j,a uxJ  -  ■  0  ' 

and  name  the  outputs  ir^ : 


The  first  of  equations  (2)  are  instrumented  and  the  outputs 

named  : 


Note  that  it  is  necessary  to  generate  In  /j. 

Since  the  outputs  of  (28)  and  (29)  are  in  the  inputs  of  (29) 
and  (28),  respectively,  it  is  sufficient  to  multiply  by  the 
appropriate  coefficients  and  interconnect;  each  is  the  feed¬ 
back  for  the  other.  The  parameters  b^  and  c,  are  available 
as  constants  or  as  functions  of  time. 

As  an  example  of  the  above  discussion,  the  N^  (nitrogen) 
distribution  lends  itself  well  since  it  is  interconnected 


with  the  other  species  only  through  ,  which  we  may  consider 
constant  for  the  present. 

With  reference  to  Model  I,  the  following  diagram  was 
•ons true  ted : 


(30) 


The  generators  for  In  are  not  shown.  The  constants  and 
satisfy 


°7  ‘  °19  "  0  ’ 

for  the  reasons  given  in  the  text;  automatically  takes  the 
proper  sign;  ♦  In (xy/  x^)  is  a  constant  since  the  concentra¬ 
tion  of  in  alveolar  sac  air  is  a  postulated  constant  and 
c  ^  is  the  tabulated  rolubllity  constant  for  N^.  This  closed 
loop  immediately  computes  the  equilibrium  Xj  and  .  At 
equilibrium,  the  rum  of  the  inputs  to  each  amplifier  is  zero, 
as  required  of  amplifiers  in  the  high-gain  mode,  and  hence 
equation  (25)  i"  satisfied.  Equations  Xj  }  0  may  be  invoked,- 


41 
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as  required,  by  means  of  a  limiter  around  the  x^  amplifier  . 
Actually  >  0,and  the  limiters  function  only  during  the 
turn-on  phase  of  the  computer.  Pinally, 


[*J 

In  — —  <  c,  ■  In  K  , 

o7]  ; 

where  the  square  brackets  mean  "concentration"  and  the  mas.- 
action  law  constant  ir>  satisfied. 

Computation  of  the  In  Xj  is  not  a  straightforward  pro¬ 
cedure.  It  1 s  done  here  by  a  linear  approximation,  but 
obviously  more  sophisticated  technique^  ^ould  be  used.  AJ.o, 
it  is  sometimes  necessary  to  attenuate  a  natural  high- 
frequency  phenomenon  that  occurs  in  high-gain  loops.  Very 
small  capacitors  may  be  used  without  affecting  the  computa¬ 
tion. 

The  basic  system  requires  Just  n  ♦  m  amplifiers;  but 

because  of  scale-factor  problems,  trace  elements  must  be 

computed  from  the  ma3s-action  constants  using  multipliers  as 

( 5) 

in  White,  et  al.,  or  by  some  other  means. 


THE  CONTINUOUS  ANALOQ  SOLUTION 

The  simplest  continuous  so’ution  for  the  equilibrium 
problem,  when  0  ■  0(Y,Tr,t),  can  be  illustrated  by  means  of 
the  example  given  above. 
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Let 

b?  =  Mt),  A  ♦  B  sin  wt.  (}2) 

If  w  is  a  moderate  frequency  —  say,  the  breathing 
rate  —  then  the  computation  ra*;e  of  the  machine  is  o  high 
that  for  any  t  the  Instantaneous  equilibrium  conditions  will 
obtain.  The  function  Q  will  then  oscillate  with  angular 
frequency  w,  a  continuous  function  of  b^(t).  Experiments  of 
this  sort  have  been  performed. 

Of  course,  each  of  the  parameters  is  available  for  manual 
investigation  of  its  effect  on  the  system.  However,  the  more 
interesting  case?  arise  when  one  attempts  to  define  the  delay 
functions,  diffusion  equations,  temperature  dependencies,  and 
such  complexities  as  may  actually  occur.  As  an  example, 
there  is  the  question  of  the  equilibrium  trajectory  when  the 
pressure  is  gradually  reduced  and  then  the  atmosphere  is 
suddenly  switched  to  pure  oxygen.  Such  a  complex  question  has 
not  been  investigated  but  seems  perfectly  feasible  of  solution. 

Methods  for  continuous  optimization  of  a  complex  system 
have  had  a  relatively  short  history  in  the  literature.  In 
general,  this  problem  is  related  to  the  synthesis  of  optimal 
control  systems  or,  better,  the  optimal  synthesis  of  optimal 
control  systems.  In  such  a  problem,  a  complex  system  under¬ 
going  transformation  is  automatically  to  be  optimised  with 
respect  to  a  quality  of  the  output.  Ordinarily,  then,  a 
feedback  link  is  required  to  alter  one  or  more  inputs,  and  a 
criterion  is  required  to  measure  the  performance  of  the  system 


with  respect  to  a  (an  optimal)  standard.  In  particular,  see 

4 

Ref.  21,  ' 

The  present  problem  is  somewhat  simpler  than  the  complete 
control  system,  but  involves  more  complex  aspects  as  the 
physiological  model  increases  in  size.  For  the  present,  we 
continuously  (in  time)  compute  the  equilibrium  values  of  a 
complex  chemical  system  as  the  inputs  vary  at  random.  Later 
on,  the  inputs  will  be  outputs  computed  from  other  physio¬ 
logical  subsystems,  and  the  whole  will  be  optimized  according 


to  viable  criteria. 
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